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ABSTRACT 


Fortran  computer  programs  hove  been  written  that  will  determine  the  geodesic  paths 
on  an  arbitrary  surface  of  revolution.  The  programs  can  also  determine  the  number 
of  circuits  of  the  geodesics  necessary  to  produce  a  wrap  of  a  specified  thickness. 
This  thickness  can  be  for  one  geodesic  or  be  the  cumulative  buildup  of  many  geo¬ 
desics.  Once  the  geodesic  paths  are  determined,  thickness  profi le  and  helix  angle 
plots  are  produced.  In  addition,  routines  are  available  for  plotting  the  geodesic 
paths  on  the  developed  surface  giving  a  two-dimensional  picture  of  the  paths  on  the 
surface . 
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SUMMARY 


Compul-er  routines  have  been  developed  for  computing  a  geodesic  path  on  an  arbi¬ 
trary  surface  of  revolution.  This  computation  is  accomplished  by  approximating  the 
surface  with  a  series  of  conical  and  cylindrical  sections  (approximating  the  contour 
of  the  surface  by  straight-line  segments)  and  determining  the  geodesic  path  on  each 
section.  It  was  shown  by  S.  P.  Gold(l)  that  the  geodesic  on  the  approximated  sur¬ 
face  converges  to  the  geodesic  on  the  actual  surface  as  the  surface  approximation 
converges. 


In  cylindrical  coordinates,  r  =  kz  +  b  on  each  of  the  sections,  and  a  geodesic  path 
can  be  written  in  terms  of  9  as  a  function  of  z .  If  the  contour  of  the  surface  is  ap¬ 
proximated  by  straight  lines  joining  the  points  (r  ,  z  )  and  the  initial  conditions  of 
the  geodesic  are  that  it  pass  through  the  point  (rQ,  zQ/  0o)  helix  angle  aQ/ 

theta  (mandrel)  rotation  (R^,)  for  one  circuit  is  found  to  be; 


R 

c 


2 


L 

L 

n  =  J 


A0  , 
n 


where  J  and  L  are  the  sections  in  which  the  geodesic  turns  around  and: 

A9_  = 

■1 


'n  =  ®  (^") 


sec  (  r  ,  ,  /  c  )  -  sec  j  r  /  c  ) 

1"+’  /  V"  /J 


on  conical  sections 


0)  /  n  J,  L, 


A9  = 

n 


+  1 


c  Jon  cylindrical  section 


ion  =  0^  , 


-1 


sec  I  ^  ^  /c)|  on  section  J, 


0  -  sec 


-1 


)] 


on  section  L, 


where: 


k=(r  ,i-rl/(z  ,,-z)  section  slope,  and 
n\n+l  n/  \n+l  n  I  ' 


c  =  rQ  sin  Q!q  . 


7 


The  geodesic  turns  around  at  the  points  where  the  surface  radius  equals  c;  that  is. 


r  .  =  c 

mm 


=  rQsin^o  . 

Thevalue  (in  radians)  when  divided  by  2  ir  is  thenumberof  revolutions  per  circuit 
for  the  specified  geodesic.  If  is  written  as  a  fraction,  A/B,  where  Aand  Bhave 
no  common  factors,  then  A  is  the  number  of  revolutions  per  pattern  and  B  the  number 
of  circuits  per  pattern.  A  routine  is  included  to  find,  if  desired,  a  new  value  of  a- 
which  will  produce  a  geodesic  with  a  specified  number  of  revolutions  per  circuit. 

Computer  routines  have  also  been  written  to  compute  and  plot  two  factors  used  in 
the  stress  analysis  of  filament-wound  structures,  helix  angle,  and  thickness  of  wrap. 
In  addition,  the  routines  will  determine  the  number  of  circuits  necessary  to  produce 
a  specified  thickness  at  a  point.  This  specified  thickness  can  either  be  from  one 
geodesic  or  be  the  cumulative  buildup  of  many  geodesics. 

One  of  the  useful  by-products  of  approximating  a  surface  by  a  series  of  conical  and 
cylindrical  sections  is  that  cones  and  cylinders  are  developable;  that  is,  if  sliced, 
they  can  be  laid  out  flat  in  a  plane  (see  Appendix  A).  Geodesics  on  a  cone  or  cy¬ 
linder  become  straight  lines  on  the  developed  surface.  Routines,  both  Fortran  and 
APT,  have  been  written  to  draw  the  developed  surface  and  to  plot  geodesic  paths 
on  this  developed  surface.  This  developed  surface  plot  has  been  useful  in  determin¬ 
ing  certain  characteristics  such  as  thickness  of  wrap  and  number  and  location  of 
crossovers.  The  developed  surface  plot  can  also  be  used  to  set  up  a  winding  machine 
by  cutting  out  the  plot  and  pasting  it  on  the  mandrel  to  be  wrapped. 

These  routines  are  useful  to  engineers  in  designing  wrap  patterns  for  filament-wound 
structures.  They  are  also  the  basis  for  routines  used  in  locating  the  path  of  a  fila¬ 
ment  feed  eye  of  a  numerically  controlled  filament-winding  machine.  (2) 


INTRODUCTION 


Combining  high-strength  filaments  with  resins  in  a  composite  structure  has  led  to 
structural  elements  and  parts  which  have  exceedingly  high  strength-to-weight 
ratios.  New  materials,  which  lend  themselves  to  filament  windings,  are  being 
rapidly  developed  and  new  applications  of  composite  structures  are  appearing.  Pro¬ 
ducts  currently  made  by  filament-winding  techniques  range  from  light-weight  fish¬ 
ing  rods  to  large  railway  tank  cars. 

As  the  applications  of  filament  winding  increase,  so  does  the  need  for  abetter  under¬ 
standing  and  definition  of  wrapping  patterns.  One  large  class  of  filament-winding 
applications  Involves  shapes  which  are  surfaces  of  revolution.  Since  a  geodesic  path 
on  any  surface  is  a  stable  path,  geodesics  are  often  chosen  as  the  desired  filament 
paths.  Therefore,  this  investigation  was  made  by  Y-12  Plant  personnel  to  determine 
geodesic  paths  on  an  arbitrary  surface  of  revolution  and  to  compute  fiber  helix  angle 
and  thickness  buildup  which  would  result  from  wrapping  these  patterns.  The  project 
was  sponsored  by  Sandia  Livermore  and  carried  out  under  Purchase  Order  ASB 
92-1849. 
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DISCUSSION  OF  THE  STUDY 


FILAMENT  PATH  ON  A  SURFACE  OF  REVOLUTION 

Since  a  geodesic  on  a  surface  is  a  stable  path,(3)  q  filament  laid  along  a  geodesic 
will  have  no  tendency  to  side  slip.  For  this  reason,  geodesics  are  often  chosen  for 
the  desired  filament  paths.  However,  the  equations  for  geodesics  on  surfaces  other 
than  simple  surfaces  such  as  spheres,  cones,  and  cylinders  are  not  easily  determined. 
Therefore,  a  method  of  approximating  a  geodesicon  an  arbitrary  surface  of  revolution 
is  undertaken. 

To  determine  a  geodesic  on  an  arbitrary  surface  of  revolution,  first  approximate  the 
contour  of  the  surface  by  a  series  of  short,  straight-line  segments.  When  rotated 
about  the  axis  of  revolution,  these  line  segments  generate  a  series  of  conical  and 
cylindrical  sections  that  approximate  the  surface  of  revolution.  Then,  by  using  the 
equations  for  a  geodesic  on  cones  and  cylinders  and  by  determining  the  criteria  for 
crossing  from  one  section  to  another,  a  geodesic  can  be  computed  for  the  arbitrary 
surface  of  revolution. 


GEODESIC  ON  A  CONE 

The  problem  associated  with  surfaces  of  revolution  can  be  simplified  by  using  cylin¬ 
drical  coordinates  (r,  z,  0).  On  a  surface  of  revolution,  r  is  a  function  of  z  and  a 
point  or  curve  on  the  surface  can  be  defined  in  terms  of  two  variables,  z  and0. 

To  determine  a  geodesic  on  a  cone,  the  property  of  the  geodesic  that  is  utilized  is 
that  between  any  two  points  on  a  surface,  the  path  of  minimum  arc  length  is  a  geo¬ 
desic.  Therefore,  to  determine  a  geodesic  between  two  points  (Figure  1),  it  is  ne¬ 
cessary  to  find  the  curve  which  minimizes  the  following  integral  (arc  length): 


j  +  (dr/dz)^  +  r^  (d9/dz)^  dz. 


A  necessary  condition(^/5)  for  the  integral  to  be  a  minimum  is: 


(1) 


_ r^  (de/dz) 

+  (dr/dz)^  +  r^  (d9/dz)^ 


d 

dz 


=  0  ,  or: 


(2) 
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Figure  1.  GEODESIC  PATH  ON  A 
CONE. 
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If  the  geodesic  passes  through  the  point  (rQ,  zq,  0q)  at  helix  angle  (a  common  way 
of  specifying  the  initial  conditions  for  a  geodesic),  the  constants  c  and  d  are  found 
(Appendix  B)  to  be; 


c  =  r^sino^Q,  and 


(Vl  +  kVk)  [o  - 


sec  (rQ/c) 


+  Bq  . 


d  =  ^>1  +  k  /k 
Thus,  the  equation  for  a  geodesic  on  a  cone  is: 

9  (z)  =  ^^1  +  k^/k^  |sec  ^  [r(z)/rQ  sin  a^]  -  (tt/2-q;q)|  +0^  . 

In  the  special  case  of  a  cylinder,  where  r=  r^,  the  differential  equation  is: 

Tq  (d0/dz)  /  ^1  +  r^  (d0/dz)^  =  c  ,  or: 


(8) 


(9) 


d9/dz  =  c  /  (r^  ^  . 

The  equation  for  the  geodesic  on  a  cylinder  becomes: 

9fe)  =  c  V^o  - 

=  (z  -  Zq)  (1  /r^)  tan  0!q  +  0q  ,  (10) 

where,  again: 

c  =  r^sina^  . 

SECTION-CROSSING  CRITERIA 

A  surface  composed  of  two  cones  is  shown  in  Figure  2. 

The  equation  for  the  surface  is: 

>■(2)  =  k^  (z  -  z^  +  r^  for  Zq  s  z  ^  z^ ,  and 
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Figure  2.  GEODESIC  ON  A  SU R FACE 
COMPOSED  OF  TWO  CONES. 

r  (z)  =  l<2  (z  -  )  +  r.j  for  <  z  <  Z2  . 

From  Equation  3  it  is  seen  that  the  geodesic  between  (r^,  Zq,  0q)  and  (r^,  z^,6.|) 
satisfies:  « 

r^  (de/dz) 

C  1  ! 


and  satisfies: 

r2  (de/dz) 


between  (r.|,  z^,  0  and  {r^,  7.^ 


(4  5) 

The  Weierstrass-Erdmann  Corner  Condition  '  is  used  to  determine  the  necessary 
crossing  condition  for  maintaining  a  geodesic  on  the  composite  surface;  that  is: 


fl  +  r^  (de/dz  +  (dr/dz)^ 


(d0/dz)  +  (dr/dz)^ 
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(It  is  shown  in  Appendix  B  that  the  condition  =  c,  implies  that  the  helix  angle  is 


continuous  at  z  =  .) 


Then  the  equation  for  the  geodesic  is: 


0(z)  =  ^  [r  (z)/c]  -  (77/2  Q!q)|  +9^ 


for  Zq  <  z  <  7.y  and 


0(z)  -  ^1  +  /k^  I  sec  ^  [r(z)/c]-sec  \r^/c)|+e^ 


for  z^  <  z  <  Z2,  (12) 


where: 


e,  =9U,)  =  ^i  tk^Vk,) 


and; 


c  =  roSmao- 


-1 


sec  (r,/ c)  -  (tt/2  -  a  ) 

I  o 


It  was  shown  (Equation  6)  that  a  geodesic  on  any  cone,  r(z)  =  k^z  +  b  , 
satisfies  the  differential  equation: 


de/dz  =  c^  Vi  +  ) 


If  the  geodesic  is  a  continuation  of  a  geodesic  which  passed  through  the  point  (rQ, 
Zq,  Qq)  helix  angle  the  constant  of  integration,  c^,  is  (see  Equation  11  ): 


c  =  c 
n 


=  roSin 


By  rewriting  the  differential  equation  as: 


dz/  d9  =  r  ^r^  -  c^  /  ^c  ^1  + 
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It  is  immediately  seen  that: 


dz/  de 


Thus,  the  turnaround  point  of  the  geodesic  is  that  location  where  the  radius  of  the 
surface  equals  c  (ie,  equals  rQ  sin  q;q);  that  is,  the  radius  at  the  turnaround  is  deter¬ 
mined  by  rQ  and  QfQ  (radius  and  helix  angle  at  the  initial  point)  and  is  independent 
of  the  shape  of  the  surface. 


DETERMINING  A  CIRCUIT  OF  THE  GEODESIC 


Let  the  contour  of  the  surface  be  defined  by  a  series  of  straight-line  segments  join¬ 
ing  the  points  (r  ,  z^),  n  =  1 ,  .  .  M  (Figure  3).  For  each  segment,  define  the 
parameter  k^,  as  follows: 


1,  2,  .  .  .,  M  -  1  . 


If  the  initial  conditions  for  specifying  the  geodesic  are  that  the  geodesic  must  pass 
through  a  point  Pq,  whose  radius  is  rg  at  helix  angle  Wq,,  then  the  constant  of  in¬ 
tegration,  c,  is: 


c 


sin  0! 


0  • 
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Note  here  that  c  must  be  such  that: 


c  > 


jmox  r,, 


/ 


otherwise,  the  geodesic  would  continue  beyond  the  defined  portion  of  the  surface. 
Now,  to  determine  the  sections  In  which  turnaround  occurs.  It  is  necessary  to  find 
J  and  L  such  that: 


,  kj  >  0, 

and 

''L  +  l 

/  k|^  <  0. 

Section  J  will  be  called  the  lower  turnaround  section  and  L  the  upper  turnaround 
section.  When  n J,L: 


Ae 


sec''  ,  /  cj-sec'' 


If  k  0 


(conical  section) 


("n  +  l  -"n)  =/('■„  V^7^)  = 


(cylindrical  section) 


(13) 


When  n  =  J,  L, 


A9j  =  (V'  +  /  ''j)  ('■j+  1  /<=)-0 


,  and 


(14) 


ABl  =  (Vl+k2/\)  0-sec-' 


(15) 
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The  rotation  during  one  circuit,  R^,  becomes: 


R  =2 
c 

In  order  for  the  geodesic  to  return  to  its  starting  point  (ie,  complete  one  pattern), 
R^  (In  revolutions)  must  be  a  rational  number,  say  R^  =  A/B.  (In  practice,  R^ 
will  always  be  rational  since  it  Is  a  computed  value.)  Then,  after  B  circuits,  the 
mandrel  will  have  completed  A  revolutions  and  the  geodesic  will  have  returned  to 
its  starting  point.  If  A  and  B  have  common  factors,  the  geodesic  will  return  to  its 
starting  point  after  fewer  circuits.  Thus,  to  determine  when  the  path  starts  re¬ 
peating,  it  is  necessary  to  reduce  A/B  to  a  fraction  which  has  no  common  factors. 
Once  this  Is  done,  A  becomes  the  number  of  revolutions  per  pattern  and  B  the  number 
of  circuits  per  pattern. 

Often  the  initial  helix  angle,  Q;q,  is  only  an  estimate  of  the  desired  helix  angle 
at  Pq.  It  may  be  more  desirable  to  have  a  helix  angle  approximately  equal  at 
Pq,  but  which  will  produce  a  wrap  having  a  predetermined  number  of  circuits  per 
pattern.  This  is  the  case  when  complete  coverage  is  desired  at  a  given  parallel  or 
where  a  certain  thickness  is  wanted  at  a  parallel.  In  Appendix  B,  an  iterative 
scheme  for  choosing  a  new  value  foi"  0(0  is  derived  to  achieve  the  number  of  circuits 
per  pattern. 

GEODESIC  ON  A  DEVELOPED  SURFACE 

One  of  the  useful  by-products  of  approximating  a  surface  of  revolution  by  a  series 
of  conical  and  cylindrical  sections  is  that  cones  and  cylinders  are  developable. 
That  is.  If  sliced,  they  can  be  laid  out  flat  in  a  plane.  To  further  simplify  matters, 
geodesics  on  a  cone  or  cylinder  become  straight  lines  on  the  developed  surface. 
Thus,  a  two-dimensional  picture  of  a  geodesic  on  the  surface  can  be  drawn. 

Drawing  a  geodesic  on  a  developed  surface  has  been  helpful  in  determining  certain 
characteristics  of  a  geodesic  such  as  the  thickness  of  the  wrap  and  the  number  and 
location  of  the  crossovers.  The  developed  surface  plot  could  also  be  used  in  setting 
up  a  winding  machine  by  cutting  out  the  plot  and  pasting  it  on  the  mandrel  to  be 
wrapped. 

Computer  routines  have  been  written  to  compute  a  geodesic  on  a  surface  for  given 
Initial  conditions,  to  develop  the  surface,  and  to  plot  the  geodesic  on  the  developed 
surface.  As  an  example  of  this  plot,  geodesics  were  computed  for  the  surface  shown 
in  Figure  4.  The  initial  helix  angles  were  adjusted  (by  the  scheme  discussed  in  Ap¬ 
pendix  B)  so  that  the  geodesic  had  11  circuits  per  pattern  (thus  returning  to  its  start¬ 
ing  point  after  1 1  circuits).  Figure  5  shows  a  single  geodesic  on  the  surface.  Figure 
6  is  the  two-dimensional  picture  of  the  geodesic  on  the  developed  surface,  and 
Figures  7  and  8  show  the  combined  pattern  of  four  geodesics  on  the  surface. 


L 

E 

n  =  J 


A6 


(16) 
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Figure  4.  CONTOUR  OF  A  SURFACE  COMPOSED  OF  TWO 
CONICAL  SECTIONS  AND  A  CYLINDRICAL  SECTION. 


B 

Figure  5.  SURFACE  WITH  A  SINGLE 
GEODESIC. 

GEODESIC  CHARACTERISTICS 

Two  parameters  used  In  the  stress  analysis  of  a  filament-wound  structure  are  the  helix 
angle  and  thickness  of  the  wrap  at  various  parallels.  The  helix  angle  can  be  de¬ 
termined  directly  from  the  relationships  (Appendix  B): 


tan  a,  = 


OL  = 


where: 


In  determining  the  thickness  of  wrap  at  a  given  parallel,  it  is  assumed  that  the 
center  of  the  band  follows  the  geodesic  path.  The  approach  used  is  to  determine, 
at  the  desired  parallel,  the  percentage  of  the  circumference  covered  by  a  circuit 
of  the  geodesic.  If  the  circuits  are  uniformly  spaced  around  the  part,  then  the 
computed  percentage  of  coverage  can  be  used  to  determine  the  average  thickness 
at  that  parallel;  that  is: 
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Figure  6.  DEVELOPED  SURFACE  WITH  A  SINGLE  GEODESIC. 


Average  Thickness  at  a  Parallel 

=  (coverage/circuit)(number  of  circuits)(band  thickness).  (18) 

For  the  derivation  of  the  equations  for  coverage  at  a  parallel,  see  Appendix  C. 

It  should  be  noted  here  that  the  value  computed  for  thickness  is  actually  the  amount 
of  glass  at  the  parallel.  It  does  not  take  into  account  the  matrix  material  present 
or  the  thickness  resulting  from  voids  and  bridging  of  the  fibers.  Therefore,  this 
figure  should  be  modified  by  some  factor  determined  by  the  percent  glass  of  the 
wrap. 


Figure  7.  SURFACE  WITH  FOUR  GEO¬ 
DESICS. 


Figure  8.  DEVELOPED  SURFACE  WITH  FOUR  GEODESICS. 


For  a  given  geodesic.  Equation  18  con  be  used  to  determine  the  number  of  circuits 
necessary  to  build  up  a  desired  thickness  at  a  parallel: 


20 


Number  of  Circuits  = 


(desired  thickness  at  the  parallel) _ 

(coverage/circuit  at  paral lel)(band  thickness) 


However,  knowing  the  number  of  circuits  to  be  wrapped  does  not  fully  describe  the 
wrap  pattern.  It  may  be  desirable  to  have  these  circuits  uniformly  spaced  around  the 
part.  This  possibility  brings  up  an  interesting  question:  Of  how  many  patterns  and 
circuits  per  pattern  should  the  wrap  consist?  In  trying  to  answer  this  question,  two 
approaches  are  taken.  They  appear  as  options  in  the  computer  program  (subroutine 
NOCIRC,  described  in  Appendix  D). 


Option  1  -  When  the  surface  to  be  wrapped  is  primarily  a  cylinder,  it  may  be 
desirable  to  have  the  circuits  spaced  around  the  part  so  that  after  one  pattern, 
the  cylindrical  portion  is  completely  covered.  Here,  the  number  of  circuits 
per  pattern  is  chosen  to  give  complete  coverage  at  a  parallel  with  no  over¬ 
lapping  of  fibers  going  in  the  same  direction.  The  number  of  patterns  neces¬ 
sary  to  build  up  the  desired  thickness  is  then  determined. 


Option  2  -  When  wrapping  a  general  surface  of  revolution,  complete  coverage 
at  one  parallel  would  produce  overlapping  fibers  or  less  than  complete  cover¬ 
age  at  all  other  parallels.  Therefore,  it  is  felt  that  the  idea  of  complete 
coverage  at  a  parallel  has  less  meaning  here.  Also,  in  wrapping  a  general 
shape,  it  may  be  desirable  to  apply  many  different  geodesics,  building  up  a 
thin  layer  with  each  to  achieve  an  overall  wrap  of  a  given  thickness.  The 
different  geodesics  could  be  chosen  to  produce  this  wrap.  Thus,  with  this 
option,  the  number  of  circuits  per  pattern  is  chosen  to  equal  the  total  number 
of  circuits  to  be  wrapped  for  the  geodesic.  Hence,  after  one  pattern,  the  de¬ 
sired  thickness  for  that  geodesic  is  obtained. 

To  achieve  a  desired  thickness  at  a  parallel,  the  number  of  circuits  per  pattern, 
and  number  of  patterns  are  determined  by  use  of  one  of  the  two  options.  The  de¬ 
sired  thickness  could  be  for  this  particular  geodesic  or  the  cumulative  thickness  of 
this  and  all  prior  geodesics.  If  it  is  the  cumulative  thickness  that  is  wanted,  then 
the  thicknesses  resulting  from  the  previous  geodesics  are  computed  and  subtracted 
from  the  thickness  specified.  This  value  is  then  used  in  determining  the  desired 
number  of  circuits.  However,  the  number  of  circuits  per  pattern  of  the  geodesic 
determined  by  the  specified  initial  conditions  will  not,  in  general,  be  the  same  as 
those  needed  to  give  this  wrap.  Hence,  it  may  be  necessary  to  find  a  geodesic 
which  differs  slightly  from  the  initially  specified  one,  but  which  has  the  needed 
number  of  circuits  per  pattern. 


The  procedure  for  finding  the  new  path  is  as  follows:  If  A/B  is  the  computed  re¬ 
volutions  per  circuit  of  the  specified  geodesic  and  NB  the  desired  circuits  per 
pattern,  an  integer  NA  is  found  so  that  NA/NB  is  as  close  as  possible  to  A/B. 
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If  NA  and  NB  have  common  factors,  NA  and/or  NB  are  altered  so  that  there  are 
no  common  factors.  Then  NA  becomes  the  number  of  revolutions  per  pattern  and 
NB  the  circuits  per  pattern.  A  new  geodesic  having  NA/NB  revolutions  per  circuit 
can  then  be  found  (by  the  scheme  described  In  Appendix  B)  or  the  rotation  of  the 
computed  geodesic  can  be  distorted  to  achieve  the  desired  revolutions  per  circuit. 


Computer  routines  have  been  written  for  plotting  these  geodesic  characteristics 
(helix  angle  and  thickness).  For  plotting  purposes,  distance  along  the  contour  of 
the  surface,  S,  was  chosen  as  the  reference  (see  Figure  4).  For  consistency  and 
ease  in  plotting,  all  of  the  quantities  are  normalized  before  plotting. 


Plots  were  made  for  the  surface  and  geodesics  shown  in  Figure  7.  Values  were  com- 
puted  for  a  0. 6-Inch-wIde  band,  0.01  Inch  in  thickness.  The  first  plot.  Figure  9, 
relates  R  and  Z  to  the  reference  S;  Figure  10  is  a  plot  of  the  helix  angles  for  the 
four  geodesics.  Figure  11  is  the  thickness  plot  for  one  geodesic  (the  geodesic  shown 
in  Figure  5),  Figure  12  shows  the  thickness  resulting  from  the  four  geodesics,  and 
Figure  13  is  a  scale  drawing  of  the  contour  after  the  wrap. 


SMAX  =  16.1933 
RMAX  =  4.0000 
2MAX  =  14.7000 


Figure  9.  PLOT  R EL ATIN G  R  AND  Z  (NORMALIZ ED)  TO  RE  PER EN  CE  S. 


S  /  SMAX 
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COMPUTER  PROGRAMS 

The  programs  for  computing  a  geodesic  and  plotting  its  characteristics  are  written  in 
Fortran  II.  There  are,  in  addition,  four  APT  macros  available  for  computing  a  geo¬ 
desic  and  plotting  it  on  the  developed  surface. 

Fortran  Program 


The  Fortran  program  consists  of  two  main  programs  and  17  subroutines.  In  addition, 
the  plotting  routines  utilize  several  subroutines  for  the  Gerber  Scientific  Plotter.  (6) 
With  slight  modification,  the  Gerber  subroutines  could  be  used  with  other  plotting 
machines. 

The  geodesic  subroutines  are  called  by  one  of  the  main  programs.  Main  program 
DESIGN  is  utilized  when  computing  and  plotting  geodesic  characteristics;  main 
program  DEVPLT  is  used  for  plotting  geodesics  on  a  developed  surface.  Flow  sheets 
of  the  main  program  and  deck  arrangements  for  the  two  operations  are  shown  in 
Appendix  D.  Also  given  in  Appendix  D  are  input  details  and  a  listing  of  the  com¬ 
puter  program. 

APT  Program 


The  APT  program  represents  the  initial  efforts  on  this  project.  Due  to  the  limited 
amount  of  storage  available  in  APT,  this  approach  was  abandoned  and  the  Fortran 
program  undertaken.  Therefore,  the  APT  program,  consisting  of  four  macros,  is 
limited  to  computing  a  geodesic  and  plotting  it  on  the  developed  surface.  These 
macros  are  briefly  described  in  Appendix  D. 


COMPARISON  OF  A  TRUE  GEODESIC  WITH  A  GEODESIC  COMPUTED  BY  THE 
APPROXIMATION  TECHNIQUE 

The  technique  described  in  this  report  is  the  computation  of  a  geodesic  for  a  surface 
which  is,  in  effect,  an  approximation  of  some  other  surface.  A  logical  question  to 
be  raised  is  just  how  good  does  this  computed  path  conform  to  a  geodesic  on  the 
original  surface?  S.  P.  Gold  proves  that  the  path  on  the  approximated  surface  con¬ 
verges  to  the  geodesic  on  the  true  surface  as  the  surface  approximation  converges. 

As  an  example  of  how  well  the  approximation  technique  works,  geodesics  on  a 
sphere  were  compared  to  those  computed  by  the  approximation  technique.  A  fila¬ 
ment  will  be  on  the  mandrel  surface  even  if  a  coarse  approximation  is  used  in  cal¬ 
culating  its  path.  For  a  given  point  (r,  z,  8)  on  the  filament  path,  there  will  be, 
for  a  given  z,  no  error  in  r  (since  the  point  lies  on  the  mandrel  surface)  between  the 
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filament  path  and  the  true  geodesic  (great  circle).  The  deviation,  if  any,  will  be 
in  the  rotation,  0.  Therefore,  in  comparing  the  computed  path  with  the  great  circle, 
the  rotation  for  a  great  circle  (360°)  is  compared  with  the  rotation  as  computed. 

Geodesics  were  computed  for  six  approximations  of  the  sphere.  These  approximations 
ranged  from  18  conical  sections  (19  equally  spaced  points  on  the  sphere)  to  an  ap¬ 
proximation  involving  720  sections.  Geodesics  with  helix  angles  (at  the  equator  of 
the  sphere  )  of  10  to  85  degrees  were  determined.  The  results  are  summarizedin 
Table  1 .  It  can  be  seen  from  this  table  that  the  approximation  technique  determines 
a  path  which  closely  follows  the  true  geodesic  on  a  sphere.  The  finer  the  approxi¬ 
mation  of  that  portion  of  the  sphere  on  which  the  geodesic  travels,  the  smaller  the 
deviation  between  the  great  circle  and  the  computed  path. 


Table  1 


COMPARISON  OF  THE  ROTATION  OF  A  TRUE  GEODESIC  ON  A  SPHERE  WITH 
GEODESICS  ON  VARIOUS  APPROXIMATIONS  OF  THE  SPHERE 


Helix 

Angle 

(degrees) 

Number  of 

Conical  SecHons 
Approximating  a  Sphere 

Number  of 

Sections  Traversed 
by  the  Geodesic 

Rotation 
for  Circuit 
(degrees) 

Deviation 
per  Circuit 
(degrees) 

Percent 

Deviation 

40 

18 

10 

366.653 

+6. 653 

1.85 

65 

36 

10 

366.433 

+6.433 

1.79 

80 

90 

10 

366.342 

+6.342 

1.76 

85 

180 

10 

366.331 

+6.331 

1.76 

40 

36 

20 

362.244 

+2.244 

0.622 

70 

90 

20 

362.196 

+2. 196 

0.612 

80 

180 

20 

362.184 

+2.184 

0.607 

85 

360 

20 

362.185 

+2.185 

0.607 

50 

90 

40 

360.773 

+0.773 

0.215 

70 

180 

40 

360.765 

+0. 765 

0.212 

80 

360 

40 

360. 764 

+0. 764 

0.212 

85 

720 

40 

360.770 

+0. 770 

0.214 

10 

180 

80 

360.269 

+0.269 

0.075 

70 

360 

80 

360.267 

+0.267 

0.074 

80 

720 

80 

360.268 

+0.268 

0.074 

10 

180 

160 

360.071 

+0.071 

0.020 

50 

360 

160 

360.094 

+0. 094 

0.026 

70 

720 

160 

360.093 

+0. 093 

0.026 

DEFINITION  OF  TERMS 


Geodesics  -  A  path  is  called  a  geodesic  on  a  surface  if  at  each  point  of  the  path, 
the  principal  normal  coincides  with  the  normal  to  the  surface  (The  shortest  of  al I 
paths  joining  two  points  on  a  surface  is  an  arc  of  a  geodesic  .)^^) 

Meridian  -  Any  plane  which  passes  through  the  axis  of  revolution  intersects  a  surface 
of  a  revolution  along  a  pair  of  curves.  The  curves  are  called  meridians. 


Helix  Angle  -  If  P  is  a  point  of  a  geodesic  on  a  surface  of  revolution,  then  the  angle 
between  the  geodesic  and  the  meridian  at  point  P  is  the  helix  angle  at  P. 
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Pnmilfil  -  Every  plane  perpendicular  to  the  axis  of  revolution  intersects  a  surface  of 
revolution  along  a  circle,  which  is  called  a  parallel. 

Circuit  -  The  path  traced  from  a  starting  point  at  a  particular  parallel  on  a  surface 
until  the  path  crosses  the  same  parallel  going  in  the  same  direction  is  one  circuit. 

Pattern  -  The  number  of  circuits  the  path  traces  on  a  surface  in  returning  to  its  orig¬ 
inal  starting  point  is  a  pattern. 


27 


APPENDIX  A 


DEVELOPED  SURFACE 


Developing  g  Surface 

Let  the  contour  of  a  surface  be  defined  os  a  series  of  straight-line  segments  joining 
the  points  (r^,  z^),  n  =  1,  M.  Define  the  following  section  parameters: 


'n=  (n+t  ■''n)/(n 


n  =  1,  2,  . . M  -  1 


f  =  \1  +  k 
n  1  n 


n  =  1,  2,  . . M  -  1 


-  0 


X  +/z  ^-z|f 

n  I  n+  I  n  /  n 


n  =  1,  2,  . . M  -  1 


<p  =  I  k  /f  2n 
^n  I  n  '  n 


n  =  1,  2,  . . M  -  1 


R1  =  f  /  k  r 
n  n  n 


r  ifk  >0 
n  n 


r  ^  if  k  <0 
n  +  1  n 


n  =  1,  2,  . . M  -  1 


R2  =  f  /k  r 
n  n  '  n 


xc  = 
n 


X  -  R1 
n  n 

X  +  R2 
n  n 


r  ,  ,  if  k  >0 
I  n  +  I  n 

)r  if  k  <0 
[  n  n 


ifk  >0 
n 

ifk  <0 
n 


n  =  1,  2,  . . M  -  1 


n  =  1,  2,  . . M  -  1 


Utilizing  these  parameters,  the  surface  can  be  developed.  Figure  A-1  is  an  example 
of  a  surface  which  has  been  developed. 


Tranformation  of  a  Point  on  Surface  (z,  8)  to  g  Point  on  the  Developed  Surface  (x,  y) 


1 .  Find  n  such  that; 
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These  equations  reduce  to: 


(vO! 

r  - 

r  +k  (z-z  ) 

cos 

(k  /f 

)el 

\n  n/| 

n 

n  n  \  n  / 

L\" 

V  J 

4/kj  r  +1<  [z-zj  sin  iy  /  ^  )9- 

In'n/n  n\  n/  Vn^ny 


(21) 

(22) 


Drawing  a  Geodesic  on  a  Developed  Surface 


Drawing  a  geodesic  on  a  conical  section  (l<r,>0)will  be  discussed.  Since  the 
other  cases  (kn  <  0)  are  similar,  they  will  not  be  presented.  Let  the  geodesic 
enter  section  n  at  (r,  z,  9)  either  initially  or  by  transition  from  another  section. 
Define ; 


(^/'n) 

e , 

AS  =  (  k  /f  J 

AS  , 

n 

R2 

n 

if  previous  section 

was  section  n  +  1 

R1 

n 

if  previous  section 

was  section  n  -  1 

Rl 

n 

if  p,  =  R2 

1  n 

R2 

n 

if  Pi  =  Rl 

^1  n 

Case  A-1  —  9q  +  AS  when  (n  ^  J)  (Figure  A-2)  -  The  geodesic  is  the  line 

segment  between  (x^,  y^)  and  (x^,  y^),  where: 

sin 

x^  =  XC^  +  cos  +  A(3), 


Yd  =  P2  5'"  +  • 


Figure  A-2.  GEODESIC  ENTERING  AND  LEAVING  A 
CONICAL  SECTION.  (Line  on  the  Developed  Surface) 


Case  A-2  —  Sq  +  A8  ><Pn/  -*)  ^-3)  -  In  this  case,  the  geodesic 

is  represented  as  two  line  segments,  from  (x^,  y^)  to  (x^,  y^)  and  from  (x2,  to 

(xj,  y^),  as  shown  in  Figure  A-3.  With  respect  to  a  local  origin  at  (xc^,  0),  the 

line  through  (x  ,  y  )  has  the  equation: 

0  0 
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and  the  line  defining  the  end  of  the  developed  surface  satisfies  equation: 


y  cos  (p^  =  X  sin  (p^  , 


where ; 

''e'Pl  =“  ^0' 

Solving  for  the  intersection  of  the  two  lines  determines  the  point  (x^,  y^.  Then 
(^2'  ^2^  found  by: 


The  end  point  of  the  line  segment  is: 

I 

^d  ^  ^2  ®  ' 

I 

Xd  =  O2  sin  8  , 


where: 


3  =  B(,+A^-<p„ 


Translating  the  points  by  xc^  locates  the  geodesic  on  the  developed  surface. 


Case  A-3  —  8q  +  2A8  l<Pp/  when  (n  =  J)  (Figure  A-4)  -  The  geodesic  on  Section  J 

(turnaround  section)  is  again  a  line  segment  between  the  points  (Xg,  yg)  and  (xj,  y^), 
as  presented  in  Figure  A-4,  where: 


Figure  A-4.  GEODESIC  ON  A  TURNAROUND  SECTION  OF 
A  DEVELOPED  SURFACE. 


X  =  XC  +  p,  COS  jS-  , 

e  n  1  0 


Xe  =  P, 

Xj  =  xc^  +  p^  cos  +  2/^  , 

Xd  =  P]  sin  ^8q  +  2A8^  . 


The  cose,  8o  2  A)3  >  <p  ,  is  similar  to  Case  A-2  and  can  be  determined  in  a 
similar  fashion. 
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APPENDIX  B 


ADDITIONAL  GEODESIC  DERIVATIONS 


Evaluation  of  the  Constants  of  Integration 


Let  the  initial  conditions  for  specifying  a  geodesic  be  that  the  geodesic  passes 
through  point  (rQ,  zq,  0q)  at  helix  angle  aQ.  To  determine  the  constant  of  in¬ 
tegration  (c  in  Equation  3  of  the  text),  two  cases  will  be  examined. 


Case  B-1  —  Initial  Point  Lies  in  a  Cylindrical  Section  -  It  was  shown  in  Appendix 
A  that  in  transforming  a  surface  point  on  a  cylinder  (r,  z,  0)  to  a  point  on  the  de¬ 
veloped  surface  (x,  y),  the  relationship  is: 


X  =  z  +  h,  where  h  is  a  constant,  and 


Then: 


y  =  r0  . 


dy/dx  =  (dy/d0)  /  (dx/d0)  , 
=  r  (d0  / dz)  ,  or 
d0/dz  =  (1/r)  (dy/dx)  . 


But, 


dy/dx  =  tanaQ,  and 


d0/dz  -  (1/r)  tanofQ,  as  shown  in  Figure  B-1. 

The  equation  developed  for  the  geodesic  on  the  cylinder  was  determined  to  be: 


2 

r 


(d0/dz)/  /l 


+  r^  (d0/dz)^  = 


Therefore, 


c 
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Z  Y 


c  =  r  sin 

Since  the  initial  point  lies  on  the  cylinder. 


sin  at 


0  • 


Case  B-2  —  Initial  Point  in  a  Conical  Section  -  As  was  true  for  the  cylindrical 
case,  the  results  of  Appendix  A  (developing  the  surface)  will  be  utilized  here.  The 
following  relations  are  derived  from  Equation  21  and  22  of  Appendix  A: 


dy/dx  = 


d0/d2 


(^-%)  +  (dz/d9)f„  sin 

+(d^/de)f„cos  /l<„  9/f  \ 

r  (dy/dx)  cos  9/f^^  -  sin  e/f^^ 


;  (23) 


By  examining  the  initial  conditons  for  ©q  0  (no  generality  lost  here  since  a  sub¬ 
stitution  6  =  9  -  0Q  would  result  in  the  same  geodesic  shifted  by  6q),  the  initial 
point  would  appear  as  shown  in  Figure  B-2.  Again, 
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Y 


Figure  B-2.  DEVELOPED  CONICAL  SECTION. 


dy/dx  =  tan  ofQ  ,  and 


d9/di 


f  tan  a-  / 
n  0 


r  +k  (z„-z) 
n  n  I  0  r\j 


The  geodesic  on  the  conical  section  satisfies  (see  Equation  5): 


c  -  r  (d0/dz)  /  . 


The  constant,  evaluated  at  (r^,  Zq,  0),  becomes: 


■■o  ('n  V^o)  /  ' 


c  =  rgsinag 


The  constant  d  is  found  to  be: 


or 
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Helix  Angle  at  a  Parallel 

To  determine  the  helix  angle  at  a  given  parallel,  the  developed  surfacewill  again 
be  utilized.  By  examining  the  geodesic  at  0  =  0  (again  no  generality  is  lost),  it  is 
seen  that: 

=  tan  a  . 

9  =  0 
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r  sin  Ck;  =  c . 


(27) 


Equation  27  is  the  same  relationship  given  by  Clairauts'  Theorem(^)  for  a  geodesic 
on  a  surface  of  revolution. 

Since  r  is  a  continuous  function  and  the  constant  c  has  the  same  value  on  each 
section.  Equation  27  implies  that  the  helix  angle,  ot,  is  continuous  at  the  point 
of  transition  from  one  section  to  another. 


Determining  an  Initial  Helix  Angle  to  Produce  a  Geodesic  with  the  Desired  Number 
of  Revolutions  per  Circuit 


Previously,  an  equation  for  determining  the  rotation  fora  circuit  of  the  geodesic, 
Rj.,  was  derived.  This  was  given  by  Equation  16  which  is  repeated  below: 


R  =2 
c 


L 

E 

n  =  J 


Ae 


where: 


A9j=  (fj/kj 


("n+1  -"n) 


-1 


(y): 


ifk  0 

n 


ifk  =  0 
n 


c=ro  sino^Q  . 


If  the  desired  rotation  per  circuit  is  R^  (to  give  complete  coverage  or  a  desired 
thickness,  etc),  define: 


ARc  =  Rc  -  Rc 
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An  approximation  of  an  initial  helix  angle,  (Xq,  which  will  produce  a  geodesic 
having  the  desired  rotation  per  circuit  is  found  by: 


Now, 


+  Aa 


0  * 


dR  /  da  =  d  2 
c  0 

=  2  E 

n  =  J 


L 

L  Ae 


/ 


d  A6  /  d 

n  0 


The  derivatives  are  found  to  be: 


dA9  / da^  = 
n  0 


dA9j/do-Q=  Vr/  +  ,  -  ' 

dA9j^/dag=  cos  a^/  -  c  , 

^  r  COS  a  /  k  \  -l/Vr 

K  /  L 


'n  +1  + 


1 

^  n 


ifk 

n 


■■o  “o 


H' 


^n+1  ■  "^n 


,^2_c2\3/2 

*  n 


ifk  =0 
n 
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APPENDIX  C 

DERIVATION  OF  EQUATIONS  FOR  THE  THICKNESS  OF  WRAP 


Computing  the  Coverage 

As  stated  previously,  the  approach  used  in  determining  the  thickness  of  wrap  at  a 
parallel  is  to  find  the  fraction  of  the  circumference  which  is  covered  by  one  circuit 
of  the  geodesic.  In  computing  this  coverage,  it  is  assumed  that  the  center  of  the 
band  lies  along  the  geodesic.  The  developed  surface  is  utilized  in  each  of  the  four 
cases  considered  below. 

Case  C-1  Parallel  in  g  Cylindrical  Section  -  Each  circuit  of  a  geodesic  will 
cross  the  parallel  in  a  cylindrical  section  twice  as  shown  in  Figure  C-1 .  The  frac¬ 
tion  of  the  circumference  covered  by  each  circuit  is: 


COVERAGE/CIRCUIT  -  2(w/cosa)  /  2Trr 

=  w/  IT  r  cos  a , 


where: 

w  represents  the  band  width, 
r  the  radius  of  cylindrical  section,  and 
a  the  helix  angle. 


Figure  C-1.  DEVELOPED  CYLINDER. 


The  equations  for  the  band  and  the  parallel  circle  are: 

Line  1  -  (-sin  a)  x  +  (cos  a)  y  =  w/2; 

Line  2  -  (-sin  a)  x  +  (cos  a.)  y  -  -w/2;  and 
Circle  1-  (x  +  p)^  +  y2  =  p^. 

Solving  for  the  intersection  of  Circle  1  with  Lines  1  and  2  gives: 

=  (-cos  a)  (p  sin  ot  -  w/2)  +  sin  a  "^p^  -  (p  sin  ot  -  w/2)2  , 

x^  =  (cos  a/sin  O')  y^  -  w/2  sin  a  , 
y^  -  (  -  cos  O')  (  p  sin  O'  +  w/2)  +  sin  a  -  ( p  sin  a  +  w/2)^  , 

=  (cos  o/sin  O')  y2  +  w/2  sino'. 


The  angles  covered  are: 


Then, 


COVERAGE/CIRCUIT  =  2(A(Pt  +A(pJ/cp, 


=  2  (Acp^  +  A(/32)/ 


f  /k 

n  n 


k  /f 


27r , 


(Acp^  +  A(P2  )/  7T 


where : 


r  =  radius  of  surface  at  the  parallel. 
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<p=  k  /f  (2 it),  and 
n  n 

a.=  helix  angle  at  the  parallel . 


Case  C-3  —  Parallel  in  a  Conical  Section,  Band  Crosses  the  Parallel  Once  -  For 
parallels  near  the  turnaround  parallel,  the  band  will  only  cross  the  parallel  one 
time  (Line  2  of  Figure  C-3  does  not  Intersect  Circle  1).  Figure  C-4  shows  this  case. 
Again,  the  portion  of  the  developed  cone  angle  covered  is  computed. 


Figure  C-4.  DEVELOPED  CONICAL  SECTION, 
lei  Once) 


(Bond  Crosses  the  Porol- 


The  points  of  intersection  are: 

=  {-cos  a)  (p  sin  a.  -  w/2)  +  sin  a.  -  (p  sin  a  -  w/2)^  , 

x.|  =  (cos  a/ sin  a)  y^  -  w/2  sinQ!  , 

y2  =  (-cos  0!)  (p  sin  O'  -  w/2)  -  sin  a  ^p^  -  (p  sin  O!  -  w/2)^  , 
X2  =  (cos  o/sin  a)  ^2  ~  ^  * 


43 


The  angles  covered  are: 


A<P-5  -  tan 


-  tan 


=  tan  ^  ^^l)]  ' 

fan"^  if  p>  1^21' 

|p  +  X2|/y2j  + it/2  if|x2  >p. 


f  P>  X2  , 


Then, 


COVERAGE/CIRCUIT  =(A(P^  +  ACP2 


IT=/a<P,  + 

=  ^A(P,  +A<()2y|kyf^|2„, 

=  I'n  / '^n  I  (^^1  ^  ■^»’2) 


Case  C-4  —  Parallel  in  a  Conical  Section;  Parallel  Beyond  Turnaround  Parallel  - 


Since  the  band  has  a  finite  width,  parallels  beyond  the  turnaround  parallel  can  be 
covered  (turnaround  parallel  being  that  parallel  where  the  geodesic  turns  around) . 
This  case  is  illustrated  in  Figure  C-5.  The  portion  of  the  angle  covered  is  again 
computed . 
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Solving  for  the  point  of  intersection  gives: 
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APPENDIX  D 


COMPUTER  PROGRAMS 
Fortran  Program 


The  geodesic  program  consists  of  two  main  programs  and  17  subroutines.  In  addition, 
the  plotting  routines  described  here  utilize  several  subroutines  written  for  the  Gerber 
Scientific  Plotter. 

The  subroutines  are  called  by  one  of  the  main  programs.  Main  program  DESIGN  Is 
utilized  when  computing  and  plotting  geodesic  characteristics.  Main  program 
DEVPLT  Is  used  for  plotting  geodesics  on  a  developed  surface.  The  deck  arrange¬ 
ments  for  the  two  operations  are  shown  In  Figures  D-1  and  D-2.  Flow  sheets  of 
these  two  main  programs  are  shown  In  Figures  D-3  and  D-4.  The  1 7  subroutines 
are  described  In  the  sections  that  follow. 

Subrogtine  PARMET  -  This  routine  computes  various  parameters  for  the  conical  and 
cylindrical  sections  that  make  up  the  surface.  The  parameters  are  stored  and  used 
by  other  routines  called. 

Subroutine  DEVELP  -  This  routine  plots  the  developed  surface.  Certain  parameters 
computed  In  this  routine  are  utilized  by  the  subroutine  which  plots  a  geodesicon 
the  developed  surface. 

Subroutine  DELTHA  -  This  routine  determines  the  delta  theta  (mandrel  rotation)  In 
each  conical  and  cylindrical  section  and  the  total  rotation  for  one  circuit.  The 
length  of  filament  laid  down  In  each  section  Is  also  computed.  The  subroutine 
argument  Is  the  geodesic  number. 

Subroutine  NOCIRC  -  This  routine  computes  the  number  of  circuits  and  the  number  of 
patterns  necessary  to  lay  down  a  given  thickness  at  a  desired  parallel.  The  first 
argument  of  NOCIRC  Is  the  geodesic  number.  The  second  argument  Is  a  flag  speci¬ 
fying  which  of  two  approaches  should  be  used  in  determining  the  number  of  circuits 
per  pattern.  If  the  flag  equals  zero,  the  number  of  circuits  per  pattern  will  be  such 
that  one  pattern  will  give  the  desired  thickness  of  the  given  parallel.  If  the  flag  is 
one,  the  number  of  circuits  per  pattern  will  be  determined  so  as  to  produce  complete 
coverage  at  the  specified  parallel.  The  number  of  patterns  necessary  to  buildup 
the  thickness  will  then  be  computed.  Since  the  specified  Initial  helix  angle  will 
not  likely  produce  a  geodesic  having  the  desired  number  of  circuits  per  pattern, 
the  third  argument  of  NOCIRC  is  a  flag  specifying  which  of  two  options  shouldbe 
taken  in  computing  the  geodesic.  If  the  flag  Is  zero,  the  initial  helix  angle  is  ad¬ 
justed  to  produce  a  geodesic  having  the  desired  number  of  circuits  per  pattern.  If 
the  flag  is  one,  the  geodesic  is  distorted  to  obtain  the  number  of  circuits  per  pattern 
that  is  wanted. 
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Figure  D-1.  CARD-DECK  ARRANGEMENT  FOR  COMPUTING  GEODESICS  AND  PLOTTING  THEIR 
CHARACTERISTICS. 
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Figure  D-2.  CARD-DECK  ARRANGEMENT  FOR  COMPUTING  AND  PLOTTING  GEODESICS  ON  A  DE¬ 
VELOPED  SURFACE. 

Subroutine  GEOPLT  -  This  routine  plots  a  geodesic  on  the  developed  surface.  The 
arguments  of  GEOPLT  are  the  geodesic  number  and  the  number  of  circuits  to  be 
plotted. 

Subroutine  DIVSUR  -  This  routine  computes  at  surface  parallels,  the  helix  angle  at 
the  parallel  for  each  geodesic,  the  thickness  produced  by  each  geodesic,  and  the 
total  thickness  at  the  parallel.  These  values  are  written  on  magnetic  tape  for  use 
by  the  plotting  routines.  The  argument  of  DIVSUR  is  the  interval  along  the  contour 
at  which  the  above  described  values  will  be  computed. 

Subroutine  RZPLOT  -  This  is  the  routine  for  plotting  R  and  Z  versus  S  (normalized). 
The  first  two  arguments  of  the  subroutine  are  the  x  and  y  coordinates  of  the  origin 
for  the  plot.  The  third  and  fourth  arguments  are  the  lengths  of  the  x  and  yaxes, 
respectively. 
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Subroutine  ANGLPL  -  This  routine  plots  the  helix  angles  versus  S.  The  four  argu¬ 
ments  are  the  same  as  those  of  subroutine  RZPLOT. 

Subroutine  THPLOT  -  This  is  the  routine  for  plotting  the  thickness  for  a  single  geo¬ 
desic  versus  S.  This  is  a  normalized  plot  with  the  thickness normalizedwith  respect 
to  the  maximum  thickness  resulting  from  all  of  the  geodesics.  The  first  argument  of 
THPLOT  is  the  geodesic  number.  The  next  four  arguments  are  the  same  as  those  of 
RZPLOT. 
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Figure  D-4.  FLOW  SHEET  OF  THE  MAIN  PROGRAM  "DEVPLT". 


Subroutine  SUMPLT  -  This  routine  plots  the  total  thickness  (normalized)  resulting 
from  all  of  the  geodesics.  The  four  arguments  are  the  same  as  those  of  RZPLOT. 

Subroutine  CNPLOT  -  This  is  the  routine  for  plotting  the  contour  of  the  surface  as 
it  appears  after  the  wrap.  The  first  two  arguments  of  CNPLOT  are  the  origin  for  the 
plot.  The  third  argument  is  the  desired  scale  of  the  X  (or  Z)  axis;  the  fourth  is  the 
scale  of  the  Y  (or  R)  axis. 
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Subroutine  PERCOV  -  This  routine  is  called  by  subroutines  NOCIRCand  DIVSUR 
to  compute,  at  a  given  parallel,  the  coverage  per  circuit  of  the  goedesic  and  the 
helix  angle  occurring  at  the  parallel.  The  first  argument  of  PERCOV  is  the  radius 
of  the  part  at  the  parallel .  The  next  argument  is  the  number  of  the  section  in  which 
this  parallel  lies.  The  third  argument  is  the  geodesic  number.  The  routine  returns 
the  coverage  per  circuit,  which  is  the  fourth  argument,  and  the  helix  angle  at  the 
parallel,  which  is  the  fifth  argument. 

Subroutine  ADJUST  -  This  routine  is  called  by  NOCIRC  to  adjust  the  initial  helix 
angle  in  order  to  obtain  a  geodesic  having  a  predetermined  number  of  revolutions 
per  circuit.  The  first  argument  of  ADJUST  is  the  geodesic  number.  The  second 
argument  is  the  desired  number  of  revolutions  per  pattern,  and  the  third  is  the  de¬ 
sired  number  of  circuits  per  pattern.  The  fourth  argument  is  the  number  of  revolu¬ 
tions  per  circuit  of  the  geodesic  as  initially  specified.  The  fifth  argument  is  the 
maximum  difference  that  will  be  allowed  between  the  number  of  revolutions  per  cir¬ 
cuit  of  a  new  geodesic  and  the  revolution  per  circuit  desired.  The  sixth  argument 
of  ADJUST  is  a  flag  indicating  to  the  calling  program  whether  or  not  a  geodesic 
could  be  found  having  the  desired  number  of  revolutions  per  circuit. 

Subroutine  NOFACT-  This  subroutine,  called  by  NOCIRC,  checks  two  integers 
for  common  factors.  If  the  integers  are  found  to  have  common  factors,  one  (or) 
both  is  altered  to  obtain  new  integers  having  no  common  factors.  The  two  argu¬ 
ments  of  NOFACT  are  the  two  integers  involved. 

Function  IGCD  -  This  is  a  Fortran  function  for  determining  the  greatest  common 
divisor  of  two  integers.  This  function  is  called  by  subroutine  NOFACT.  The  two 
arguments  are  the  two  integers  whose  greatest  common  divisor  is  desired. 

Subroutine  SHIFT  -  This  routine  is  called  by  subroutine  DEVELP  if,  when  plotting 
the  developed  surface,  two  of  the  sections  overlap.  The  first  argument  of  SHIFT 
is  the  section  number.  The  second  argument,  computed  by  the  subroutine,  is  the 
amount  of  shift  necessary  to  prevent  the  section  from  overlapping. 

Subroutine  AXPLOT  -  This  routine,  called  by  the  various  routines  for  plotting 
geodesic  characteristics,  is  an  axis  generator.  Its  purpose  is  to  draw  and  label  the 
axes  for  a  plot.  The  first  two  arguments  are  the  x  and  y  coordinates  of  the  origin 
for  the  plot.  The  third  and  fourth  arguments  are  the  length  of  x  and  y  axes,  re¬ 
spectively.  The  fifth  and  sixth  arguments  are  the  divisions  per  inch,  on  the  x  axis 
and  y  axis,  to  be  marked.  The  seventh  and  eighth  arguments  are  the  length  each 
division  represents  (x  and  y  axes).  The  ninth  and  tenth  arguments  specify  which 
divisions  are  to  be  labeled;  ie,  if  this  number  is  one,  every  division  will  be  label¬ 
ed;  if  two,  every  other  division  will  be  labeled,  etc.  The  eleventh  argument  is  a 
flag  indicating  the  size  of  letters  to  be  used  in  labeling  the  axes.  Thetwelveth 
and  thirteenth  arguments  are  the  names  of  the  x  and  y  axes,  respectively. 
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Input  Format  for  the  Main  Program  DESIGN 


The  Input  format  for  computing  a  geodesic  and  plotting  its  characteristics  is  shown 
in  Figure  D-5.  The  first  card  of  the  input  is  a  title  card  containing  alphanumeric 
information.  The  next  card  states  the  number  of  points  to  be  used  in  defining  the 
contour  of  the  surface.  The  points  defining  the  surface  (rp,,  z  )  are  then  given  in 
order  of  increasing  z,  where  z  Is  the  axis  of  revolution  of  the  surface.  The  next 
card  contains  the  plotting  flags,  one  flag  for  each  type  plot  available.  Also  In¬ 
cluded  on  this  card  is  the  interval  along  the  surface  contour  (step  length)  at  which 
thickness  and  helix  angle  are  to  be  computed  and  plotted.  (If  all  of  the  plotting 
flags  are  zero,  thickness  and  helix  angle  will  not  be  computed  at  al I . )  The  scale 
of  the  final  contour  plot  (1.0  =  full  scale,  .5  =  half  scale,  etc  )  is  also  included 
on  this  card.  The  next  card  indicates  the  number  of  geodesics  to  be  wrapped.  Then 
finally  there  is  one  card  specifying  each  geodesic.  This  card  contains  the  initial 
point  of  the  geodesic  (rg  and/or  zq  is  needed),  the  initial  helix  angle  (this  could  be 
90°  if  the  initial  point  is  the  turnaround  point),  the  desired  thickness  at  a  specified 
point  (r  and/or  z  needed),  the  band  dimensions,  and  three  flags.  The  first  flag  in¬ 
dicates  which  option  is  to  be  used  in  determining  the  number  of  circuits  per  pattern 
(Options  1  and  2).  The  second  flag  indicates  whether  a  new  geodesic  having  the 
desired  circuits  per  pattern  is  to  be  found  (adjust)  or  if  the  computed  rotation  of  the 
geodesic  is  to  be  linearly  distorted  (distort)  to  produce  the  desired  number  of  circuits 
per  pattern.  The  final  flag  indicates  whether  the  desired  thickness  is  tobe  produced 
by  the  current  geodesic  alone  or  If  It  is  the  ci^mulative  thickness  of  this  and  all  prior 
geodesics. 


Input  Format  for  the  Main  Program  DEVPLT 


The  input  format  for  computing  geodesics  and  plotting  them  on  a  developed  surface 
is  shown  in  Figure  D-6.  The  format  for  the  surface  definition  is  identical  to  that 
used  in  main  program  DESIGN,  One  card  is  needed  to  specify  each  geodesic  to  be 
drawn.  This  card  contains  the  initial  radius  (tq),  the  Initial  helix  angle  (cVq)/  the 
advance  per  pattern,  the  number  of  circuits  per  pattern,  and  the  number  of  circuits 
to  be  drawn.  If  the  number  of  circuits  per  pattern  is  specified,  a  new  geodesic  will 
be  computed  (if  necessary)  to  obtain  one  having  the  desired  number  of  circuits  per 
pattern.  If  this  field  is  left  blank,  the  geodesic  as  specified  will  be  plotted. 

Output  of  the  Main  Program  DESIGN 

Output  of  the  main  program  DESIGN  and  its  subroutines  consists  of  the  following: 


(1)  The  alphanumeric  information  on  the  title  card; 

(2)  The  ini  tial  conditions  for  the  geodesic; 
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igure  0-5.  INPUT  FORMAT  FOR  THE  MAIN  PROGRAM  "DESIGN". 
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igure  D-6.  INPUT  FORMAT  FOR  THE  MAIN  PROGRAM  "DEVPLT". 
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(3)  Computed  data  for  the  geodesic  which  includes  the  geodesic  number,  the 
number  of  circuits  necessary  to  produce  the  desired  thickness,  the  number 
of  patterns  needed,  the  number  of  circuits  per  pattern,  a  ratio  of  integers 
which  is  the  ratio  of  revolutions  per  pattern  to  circuits  per  pattern  (le, 
the  number  of  revolutions  per  circuit),  the  computed  thickness  at  the 
specified  point  resulting  from  this  geodesic,  and  the  parallels  at  which 
turnaround  occurs; 

(4)  The  distortion  factor  (the  computed  rotation  is  multiplied  by  this  factor 
to  achieve  a  wrap  having  the  desired  number  of  revolutions  per  circuit); 

(5)  The  rotation  (delta  theta)  and  length  of  filament  to  be  laid  down  in  each 
conical  and  cylindrical  section  during  one-half  circuit;  and 

(6)  The  total  rotation  (degrees)  and  total  length  of  filament  for  one  circuit. 

Output  of  the  Main  Program  DEVPLT 

Output  of  the  main  program  DEVPLT  and  its  subroutines  consists  of  the  following: 

(1)  The  alphanumeric  information  on  the  title  card; 

(2)  The  initial  conditions  for  the  geodesic; 

(3)  The  rotation  and  filament  length  for  one  circuit; 

(4)  A  ratio  of  integers  which  is  the  ratio  of  revolution  per  pattern  to  cir¬ 
cuits  per  pattern; 

(5)  The  parallels  at  which  turnaround  occurs;  and 

(6)  The  rotation  and  filament  length  in  each  section  for  one-half  circuit. 

APT  Program 

The  APT  program  represents  the  initial  efforts  on  this  project.  Due  to  the  limited 
amount  of  storage  available  in  APT,  this  approach  was  abandoned  and  the  Fortran 
program  previously  described  was  undertaken.  Therefore,  the  APT  program  is 
limited  to  computing  a  geodesic  and  plotting  it  on  the  developed  surface.  The  in¬ 
put  to  the  program  is  a  point  definition  of  the  contour  to  be  wrapped.  The  19  re¬ 
served  words  should  be  dimensioned  at  least  as  large  as  the  number  of  points  de¬ 
fining  the  contour. 

MAC1  — This  routine,  utilizing  the  point  definition  of  the  surface  contour,  defines 
the  developed  surface  and  computes  various  section  parameters.  The  argument  of 
MAC1,M,  is  the  number  of  points  defining  the  contour. 

MAC2  —  This  is  the  routine  for  plotting  the  developed  surface.  The  argument,  M, 
is  again  the  number  of  points  defining  the  contour. 


MAC3  —  This  routine  computes,  for  an  initial  helix  angle  and  radius,  one  circuit 
of  the  geodesic.  The  initial  helix  angle  is  adjusted  (if  necessary)  to  produce  a 
geodesic  having  a  specified  number  of  circuits  per  pattern.  The  first  argument  of 
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MACS,  RO,  is  i-he  radius  of  the  initial  point.  The  second  argument,  AZERO,  is 
the  initial  helix  angle  at  the  starting  point.  The  third  argument,  PRIME,  is  the  de¬ 
sired  number  of  circuits  per  pattern  and  should  be  a  prime  number.  If  other  than  a 
prime  number  is  specified,  the  program  could  determine  a  geodesic  having  fewer 
circuits  per  pattern  than  desired.  The  fourth  argument,  M,  is  the  number  of  points 
defining  the  surface  contour.  The  final  argument,  EPS,  is  the  maximum  allowable 
difference  between  the  number  of  revolutions  per  circuit  of  the  computed  geodesic 
and  the  revolutions  per  circuit  desired.  This  value  will  normally  be  small.  How¬ 
ever,  if  the  user  does  not  want  the  initial  helix  angle  altered,  a  large  value  (say, 
EPS  =  1)  should  be  used. 

MAC4  —  This  is  the  routine  for  plotting  the  geodesic  on  the  developed  surface.  The 
first  argument,  TZERO,  is  the  starting  value  of  theta,  6g.  The  second  argument, 
J,  is  the  section  in  which  the  plot  will  originate.  The  third  argument,  NUMBER, 
is  the  number  of  circuits  to  be  drawn.  The  plot  will  begin  at  the  left  side  of  section 
J,  proceed  to  the  right,  and  terminate  at  the  right  hand  side  of  section  J-1.  If  the 
value  specified  for  NUMBER  is  the  same  as  that  specified  for  PRIME  in  MAC3,the 
geodesic  will  return  to  its  starting  point  (ie,  complete  one  pattern). 

A  limited  amount  of  program  output  appears  after  subroutine  MAC3  has  been  ex¬ 
ecuted.  The  parameters  printed  out  are: 

(1)  Pass  -  The  number  of  iterations  of  the  initial  helix  angle  to  get  desired  number 
of  circuits  per  pattern; 

(2)  Del  -  The  difference  between  the  desired  number  of  revolutions  per  circuit  and 
revolutions  per  circuit  actually  obtained; 

(3)  Rvn  -  revolutions  per  circuit  obtained  with  adjusted  initial  helix  angle; 

(4)  Integer  -  the  integral  part  of  Rvn; 

(5)  Fract  -  the  fractional  part  of  Rvn; 

(6)  Partn  -  the  fractional  part  of  desired  number  of  revolution/circuit; 

(7)  N  -  the  numerator  of  (N/PRIME)  which  results  in  the  value  of  Partn; 

(8)  Alpha  -  the  adjusted  initial  helix  angle; 

(9)  Tsum  -  the  rotation  for  one  circuit  in  degrees; 

(10)  Cons  -  the  constant  of  integration,  rg  sin  oig  (also  radius  at  turnaround); 

(11)  L  -  the  upper  turnaround  section; 
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(12)  I  -  The  lower  turnaround  section; 

(13)  Dtheta(n)  -  the  rotation  occurring  on  section  n  in  one-half  circuit  (In  degrees); 

(14)  Dbeta(n)  -  the  angle  traversed  on  the  developed  surface  of  section  n  (one-half 
circuit);  and 

(15)  Flngth(n)  -  the  length  of  filament  on  section  n  (one-half  circuit). 


Program  Listing 


The  two  main  programs  and  the  17  subroutines  previously  described  are  listed  below. 
Following  this  Fortran  list  is  a  listing  the  the  APT  macros. 
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*LABEL 


cdesign  main  design  program  for  winding  geodesics 

dimension  R ( I  odd ) »Z(  1000 ) .aK(  1000) »F( I  000 ) .X  I ( 1 000) »TiTLE(  I  2 ) , 

IRO (  I OD) t^oi  I  00 ) .ALPhaOI I  00) .CONSC I  00 ) .W( I  00 ) »D(  I  00) » THICK (  I  00) » 

2RT ( I  00)  >ZT (  I  00 ) »NC( I  00) *DTHETa( I  000) .FLNgTHI  1000) 

COMMON  M.R,Z.AK»F»XI jNOGEOD ♦RO »Z0» ALPHAO , CONS » w . d » thi ck *nc » rt  »zt , 
|SMAX,RMaX,ZMAX,THMAX,JJ,TiTLE,PI »DTHETA.FLNGTH,TSUM.FLSUM,NHIGH. 

2  NLOW  .  DISTRT  ,  ADVNCE  »  SHAFT |  ,  SHAFT2 
COMMON  AA,  EB,  CC,  DEL,  DELRHO,  NSTART 
COMMON  LLL  ,  RHOMIN,  FR  ,  TMIN 
PI  X  3.141 59265 
REWIND  8 

5  READ  INPUT  TAPE  5  ,  1000  ,  (  TITLE(K)  ,  K  #  I,  12  ) 

1000  FORMAT  (  I 2A6  ) 

READ  INPUT  TAPE  5  ,  lOlO  ,  M 
I  DIO  FORMAT  (  14  ) 

READ  INPUT  TAPE  5  »  1020  »  (  R(N)  ,  Z(N),  N  #  I  ,  M  ) 

1020  FORMAT  (  6F(2.6  ) 

READ  INPUT  TAPE  5  ,  1030  ,  ITHICK  ,  ISUM,  IANGLE,  ICON, STEP,  SCLe 
1030  FORMAT  (  413  ,  2F6.3  ) 

CALL  PARMET 

READ  INPUT  TAPE  5  *  1010  >  NOGEOD 
DO  100  I  #  I  ,  NOGEOD 

READ  INPUT  TAPE  5  »  1050  ♦  RO ( I  )  ,  ZO ( I ) , ALPHAO ( I ) , TH I CK ( I )  ,  RT(I), 

1  ZT(I),  W(I)  ,  D(I)  ,  KK  ,  LA  ,  KTHICK 

1050  format  (  2FI2.6  ,  2F6.3  ,  2FI2.6  ,  2F6.3  ,  312  ) 

write  output  tape  6  ,  1080  »  (  TITLE(K)  ,  K  «  I  ,12  ) 

1080  FORMAT  (  IHI  ,  I 2A6  ) 

write  output  tape  6  ,  1090 

1090  format  (  IMHO  INITIAL  INITIAL  HELIX  DESIRED  AT  PO 

lINT  BAND  BAND  CIRC/PAT  ADJUST  THICKNESS  / 

2  I08H  R  Z  ANGLE  THICKNESS  R  Z 

3  WIDTH  Thickness  flag  flag  flag  ) 

write  output  tape  6  ,  non  ,  ROin  ,  zoin  ,  alphaoim  ,  thickid, 

I  RT(I)  ,  ZTd)  ,  W(I)  ,  D(I)  ,  KK  ,  LA  ,  KTHICK 
I  I  DO  FORMAT  (  IHO  »  8FI0.4  ,  318  ) 

IF  (  I  -  I  )  95  ,  95  ,  10 

10  IF  (  KTHICK  )  95  ,  95  ,  15 
I  5  IF  (  ZT(  I  )  )  20  ,  20  ,  70 

20  IF  (  RT(I)  )  25  .  25  ,  50 

25  write  output  TaPE  6  ,  2000  «  I 

2000  format  (3IH0(RT,ZT)  NOT  GIVEN  FOR  GEODESIC  ,  I3,I3H  (RO,ZO)  USED  ) 
IF  (  ZO(I)  )  30  ,  30  ,  40 

30  IF  (  RO(I)  )  35  ,  35  ,  45 

35  write  output  TaPE  6  ,  2010  ♦  I 

2010  EORMAT  (  44H0STARTING  STATION  NOT  SPECIFIED  FOR  GEODESIC  ,  13  , 

I  I6H  CANNOT  COMPUTE  ) 

GO  To  IDO 
40  ZT(  I  )  tf  Z0(  I  ) 

RT( I )  «  Roll) 

GO  TO  70 
45  RT( I )  «  Roll) 

50  CONTINUE 
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DO  55  N  2  >  M 

IF  (  RT( I )  -  R(N)  )  65  »  60  »  55 

55  CONTINUE 

write  output  Tape  6  .  2020  *  i 

2020  FORMAT  (  38HOCOULD  NOT  LOCATE  (RT»ZT)  FOR  GEODESIC  ,  13  »  I 7H  THIS 
I  ONE  SKIPPED  ) 

GO  TO  IDO 
60  N  N  +  I 

65  NTH  N  -  I 

GO  TO  84 
70  CONTINUE 

DO  72  N  #  2  .  M 

IF  (  ZT( I )  -  Z(N)  )  78  >  76  .  72 

72  CONTINUE 

IF  (  RT(  I)  )  74  »  74  ,  50 

74  write  output  TAPE  6  .  2020  *  I 
GO  TO  inn 
76  N  #  N  +  I 
78  NTH  #  N  -  I 

IF  (  RT(  I  )  )  82  •  82  .  84 

82  RT(I)  ff  AK(NTH)  *  (  ZT(i)  -  Z(NTH)  )  +  R(NTH) 

84  CONTINUE 
IM I  I  -  I 
DO  88  K  «  I  .  IMI 

CALL  PeRcOV  (  RT(I)  .  NTH  .  K  ,  PERCNT  ,  HANGL  ) 

TK  ff  FLOATF  (  NC(K)  )  *  D(K)  *  PERCNT 
thick!  I  )  »  THICK!  I )  -  TK 
88  CONTINUE 

IF  !  THICK! I )  )  90.90*95 

90  NC  !  I  )  ff  0 

WRITE  OUTPUT  TAPE  6  ,  2030  .  I 

2030  format  !  I 3HOFOR  GEODESIC  ,  13  ,  87H  THICKNESS  BUILT  UP  BY  PREVIOU 
|S  LAYERS  EXCEEDS  DESIRED  THICKNESS  -  THIS  ONE  NOT  NEEDED  ) 

GO  TO  100 
95  CONTINUE 

CALL  NOCIRC  !  I .  KK  ,  LA  ) 

100  CONTINUE 

IF  !  ITHICK  +  ISUM  +  lANGLE  +  ICON  )  250  .  250  >  120 

120  IF  !  STEP  )  130  .  130  .  1 40 

130  STEP  #  .050 
140  CALL  DIVSUR  !  STEP  ) 

XO  #  0.0 
YO  ff  0.0 

IF  !  ITHICK  +  ISUM  +  lANGLE  )  210  .  210  »  142 

142  CONTINUE 
XL  5.0 
YL  i¥  5.0 

CALL  RZPLOT  !  XO,  YO.  XL  ,  YL  ) 

IF  1  lANGLE  )  160  .  160.  150 

150  CALL  ANgLPL  (  XO .  YO,  XL  ,  YL  ) 

160  CONTINUE 

IF  !  ITHICK  )  190  .  190  .  170 

170  DO  180  I  «  I  .  NOGEOD 
180  CALL  THPLOT  (  I,  XO,  YO,  XL,  YL  ) 

I  90  CONTINUE 
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IF  (  ISUM  )  2ID  »  210  »  200 

200  CALL  SUMPLT  (  XO*  YO»  XL  ,  YL  ) 

210  CONTINUE 

IF  (  ICON  )  250  ♦  250  ♦  220 

220  IF  (  SCLE  )  230  ,  230  »  240 

230  SCLE  #  .5 
240  CONTINUE 

CALL  CNPLOT  (  XO,  YO,  SCLE,  SCLE  ) 
250  CONTINUE 

END  file  8 
GO  TO  5 
END 

*LABEL 


CDEVPLT  MAIN  PROGRAM  FOR  PLOTTING  GEODESIC  ON  DEVELOPED  SURFACE 

dimension  R(  I  000 ) ,21  I  000) ♦AKI  I  COO) »F ( I  DOG  )  ,X I  (  I  000 ) *TiTLE(  12), 

I RO (  I  00)  ,20 (  100) , ALPHaCK  100) ,C0NS( I  00 ) , W (  I  DO ) , D (  I  00 ) ,Th I CK (  I  00 ) , 
2RT(  I  00)  ,2T(  I  00 ) ,NC( I  00) ,DTHETA(  1000)  ,FLN6TH( 1000) 

COMMON  M,R,Z,AK,F,XI ,N0GE0D ,R0 ,20 , ALPH AO, CONS , W , D , THI CK ,NC , RT ,2T , 
|SMAX,RMAX,ZMAX,THMAX,JJ,TITLE,PI ,DTHETA,FLNGTH,TSUM,FLSUM,NHIGH, 

2  NLOW  ,  DISTRT  ,  ADVNCE  ,  SHAFT  I  ,  SHAFT2 
COMMON  AA,  SB,  CC,  DEL,  DELRHO,  NSTART 
COMMON  LLL  ,  RHOMIN,  FR  ,  TMIN 
PI  #  3, 14159265 
REWIND  8 

ID  READ  INPUT  TAPE  5  ,  1000  ,  (  TlTLE(K)  ,  K  #  I,  12  ) 

IDOO  FORMAT  (  I2A6  ) 

READ  INPUT  TAPE  5  ,  1010  ,  M 
I  DIG  format  (  14  ) 

READ  INPUT  TAPE  5  ,  1020  ,  (  R(N)  ,  Z(N)  ,  N  #  I  ,  M  ) 

1020  FORMAT  (  6FI2.6  ) 

READ  INPUT  TAPE  5  ,  lOlO  »  NOGEOD 

CALL  PARMET 

CALL  DEVELP 

DO  200  I  »  I  ,  NOGEOD 

READ  INPUT  TAPE  5  «  1030  ,  RO(I)  ,  ALPHAO(l),  ADVDEG,  NCPERP,  NUM 
1030  FORMAT  (  3FI2.6  ,  2l6  ) 

DO  20  N  #  2  ,  M 

IF  (  RO(I)  -  R(N)  )  40  ,  30  ,  20 

20  CONTINUE 

write  output  TAPE  6  ,  2000  ,  I 

2000  FORMAT  (  52H0  COUlD  NOT  DETERMINE  STARTING  SECTION  FOR  GEODESIC  , 
I  13  ,  20H  ,  THIS  ONE  SKIPPED  ) 

GO  TO  200 
30  N  #  N  +  I 
40  NSTART  »  N  -  I 

IF  (  AK(NSTART))  50  ,  60  ,  50 
50  20(1)  ff  (  RO(I)  -  R(NSTaRT)  )  /  AK(NSTaRT)  +  2(NSTaRT) 

GO  TO  70 

60  20(1)  #  2 ( NSTART  ) 

70  CONTINUE 

CALL  DELTHA(I) 

NLOW  NLOW 
NHIGH  #  NHIGH 
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IF  (  NCPERP  )  145  »  145  •  80 

80  RVN  #  TSUM  /  360.0 
INTGR  U  RVN 

FRACT  n  RVN  -  FLOATF  (INTGR) 

NB  »  NCPERP 
ANB  #  NB 

ADVNCE  »  ADVDEG  /  (  ANB  *  360.0  ) 

N  #  I 
AN  »  I .0 

PARTn  #  1.0  /  ANB 
90  AN  #  an  +  1.0 
N  »  N  +  1 
PARTN  I  «  PARTN 
PARTn  #  an  /  ANB 

IF  (  FRACT  -  PARTn  )  I  I  0  >  1 50  ♦  1 00 
1 00  IF  (  N  -  NB  +  I  )  90  «  130  *  130 

110  IF  (  ABSF(  FRACT-PARTN)  -  ABSf(  FRACT-PARTN I )  )  130.  130  .  120 

1 20  N  «  N  -  I 

PARTN  n  PARTN I 
130  CONTINUE 
NA  #  N 

EPS  »  .000001 

NAA  »  NB  *  INTGR  +  NA 

CALL  ADJUST  (  I  ,  NAA  .  NB  .  RVN  »  EPS  .  LL  ) 

NLOW  n  NLOW 
NHIGH  #  NHIGH 

IF  (  LL  )  150  .  150  .  140 

140  WRITE  OUTPUT  TAPE  6  ,  20 1  0  .  I  »  ALPHAO(I)  ,  RO(I) 

2010  format  (  28H0  COULD  NOT  ADJUST  GEODESIC  .  13  ,  53H  .  THEREFORE  PL 

lOT  IS  FOR  GEODESIC  HAVING  HELIX  ANGLE  .  F6.3.  IIH  AT  RADIUS  »F6.3) 
145  NA  »  0 
NB  #  0 
INTGR  #  0 
ADVDEG  #  0.0 
I  50  CONTINUE 

ZLOW  (  CONS(I)  -  R(NLOW)  )  /  AK(NLOW)  +  Z(NLOW) 

ZHIGH  «  (  CONS(I)  -  R(NHIGH  )  )  /  AK(NHIGH)  +  Z(NHIGH) 

write  OUTPUT  tape  6  ,  2020  »  (  TITLE(K)  ,  X  #  I  ,  12  ) 

2020  format  (  IHI  .  I 2A6  ) 

write  output  Tape  6  .  203o 

2030  format  (  93HaGEODESlC  HELIX  AT  ADVANCE  TOTAL  FILAM 

I  ENT  integers  turnaround  STATIONS  / 

I  99H  NUMBER  ANGLE  RADIUS  PER  PAT  ROTATION  LENGTH  N  + 
3A  /  B  RADIUS  Z  LOWER  Z  UPPER  ) 

WRITE  OUTPUT  TAPE  6  ,  2040  »  I.  ALPHAO(I)  ,  RO(I),  ADVDEG  »  TSUM  , 
I  FLSUM  .  INTGR  .  NA  .  NB  .  CONS(I).  »  ZLOW  .  ZHIGH 

2040  FORMAT  (  I  HO  .  14  .  5FtD.3  .  314  .  3FI0.3  ) 

WRITE  OUTPUT  TAPE  6.2050.  (  N , DTHETA ( N ) . FLNGTH ( N ) . N#NLOW. NH I GH  ) 
2050  FORMAT  (  |H0  /  39H0  SECTION  DELTA  THETA  FILAMENT  LENGTH  / 

I  (  IH  .  14  .  2FI6.6  )  ) 

WRITE  OUTPUT  TAPE  6  .  2060  .  TSUM  .  FLSUM 
2060  format  (  8H0CIRCUIT  .  F13.6  .  FI6.6  ) 

IF  (  NUM  )  160  .  160  .  170 

160  NUM  »  I 
170  CONTINUE 


CALL  GEOPLT  (  I  ,  NuM  ) 

200  CONTINUE 
GO  TO  10 
END 

*LABEL 

cparmet  compute  section  parameters 

SUBROUTINE  PARMET 

dimension  R( I OOD ) ,Z( 1000 ) »A<( I  000) »E (  I  000 )  .X  I (  I  000  )  »TiTLE(  12), 
|RO(  100)  ,ZO(  IDO) ,ALPHaO( I  00) , CONS ( I  00 ) , W ( 1 1 0 ) ,D ( I  00 ) , Th I CK (  I  00 ) , 
2RT( 100)  ,ZT(  100) »NC( I00),DTHETa( 1000) ,FLNGTH(  1000) 

COMMON  M,R,Z,AK,F,XI ,NOgEOD,RO,ZO , ALPHAO, CONS , W , D , TH I CK ,NC , RT ,ZT 
I SMAX.RMAX ,ZMAX,ThMAX, JJ, Title, PI ,DTHETA  ,FLNGTh,TSUM,ELSUM,NHIGH, 
2  NLOW  ,  DISTRT  ,  ADVnCE  ,  SHAFT  I  ,  SHAFT2 
XI  (  I  )  it  0.0 
MM  #  M  -  I 
RMAX  f?  R(  I  ) 

DO  1200  N  »  I  ,  MM 

IF  (  ABSFI  Z(N)  -  Z(N+I))  -  .0001  )  lOlO  ,  1010  ,  1040 
1010  IF  (  R(N+|)  -  R(N)  )  1020  ,  1020  ,  1030 
1 020  AK(N)  »  -  (  I .OE  20  ) 

GO  TO  1035 

I  030  AK(N)  #  I .0  E  20 
1 035  F(N)  #  I .0  E  20 

X|(N+I)  it  X|(N)  +  ABSF(  R(N+|)  -  R(N)  ) 

GO  TO  1200 

1040  AK(N)  it  (  R(N+|)  -  R(N)  )  /  (  Z  ( N+ 1  )  -  Z(N)  ) 

IF  (  ABSFI  AK(N)  )  -  .0001  )  1050  ,  1050  ,  MOO 

I  050  AKIN)  tt  0.0 

MOO  FIN)  #  SQRTFI  1.0  +  AKIN)**2  ) 

XIIN+I)  tt  XIIN)  +  I  ZIN+I)  -  ZIN)  )  *  FIN) 

1200  RMAX  it  MAXIF  I  RIN+M  ,  RMAX  ) 

ZMAX  tt  ZIM) 

RETURN 

END 

♦LABEL 


CDEVELP 

SUBROUTINE  DEVELP 

dimension  R I  I  000 ) ,ZI 1000) ,Al^t  1000) ,Fl  I  000  )  ,X  I  I  I  000) ,TlTLEl  12), 

1  RO  I  I  00  )  ,ZOM00  )  ,  ALPhaOI  I  00)  ,  CONS  I  I  00  )  ,  W I  I  00 )  ,  D  M  00  )  ,Thi  cK  I  I  00  )  , 
2RT I  100)  ,ZTI  100 ) ,NC(  100) ,DTHETa I  I  000)  ,FLNGTHI  1000) 

dimension  R I  I  I  00 ) ,R2 I  I  00 ) ,PHI I  I  00) ,xcl  100) 

COMMON  M,R,Z,AK,F,X I .NOGEOD ,R0  ,Z0, ALPHAO  ,  CONS  ,W , D  ,THI CK  ,NC , RT  ,ZT 
|SMAX,RMAX,ZMAX,THMAX,JJ,TITLE,PI ,DTHETA,FLNGTH,TSUM,FLSUM,NHIGH, 

2  NLOw  ,  DISTRT  ,  ADVNCE  ,  SHAFT  I  ,  SHAFT2 
COMMON  AA,BB,CC,DEL,DELRH0,NSTART 
COMMON  LLL  ,  RHCMIN  ,  FR  ,  TMIN 

COMMON  Rl  ,  R2  ,  PHI  ,  XC 
SHIFT  I  #0.0 
Ml  #  M  -  I 
DO  200  N  #  I  ,  Ml 
IF  I  AKIN)  )  10  ,  40  ,  70 

10  FOK  tt  -  FIN)  /  AKIN) 

RMN)  #  RIN+M  *  FOK 


on  CO  CD  CD 


R2(N)  »  R(N)  *  FOK 
PHI (N)  «  360.0  /  FOK 
XC(N)  #  XI (N)  +  R?(N) 

IF  (  AK(N)  -  AK(N-n  -  .000001  )  30  »  30  . 

20  CALL  SHIFT  (  M  .  SHIFT2  ) 

Shift  i  tt  shift  i  +  shifT2 
30  XC(N)  #  XC(N)  +  SHIFTI 
GO  TO  200 

40  R2(N)  tt  2.0  *  PI  *  R(N) 

IF  (  AK(N-I )  )  50  .  60  »  60 

50  CALL  SHIFT  (  N  ,  SHIFT2  ) 

SHIFTI  »  SHIFTI  +  SHIFT2 
60  XC ( N)  #  X  I  (N)  +  SHIFTI 
Rl  (N)  It  XI  (N+n  +  SHIFTI 
GO  TO  200 

70  FOK  tt  FIN)  /  AKIN) 

R  I  IN)  tt  RIN)  *  FOK 
R2IN)  #  RIN+I )  *  FOK 
PHI  IN)  tt  360.0  /  FOK 
XCIN)  tt  XIIN)  -  RIIN) 

IF  I  N  -  I  )  100  *  100  »  80 

80  IF  I  AKIN)  -  AKIN-I)  -  .000001  )  100  ♦  100 

90  CALL  Shift  i  n  »  shifT2  ) 

SHIFTI  tt  shift  I  +  SHIFT2 
100  XCIN)  tt  XCIN)  +  SHIFTI 
200  CONTINUE 

TITLEI7)  tt  242565254346 

TITLEI8)  tt  472524606264 

TITLEI9)  tt  512621232560 

TITLEIIO)#  474346635360 

CALL  SETUP  I  TITLE  1 

CALL  PLOT  I  0.0  .  0.0  .  I  ♦  2  ) 

DO  300  N  id  I  .Ml 
IF  I  AKIN)  )  220  »  250  »  260 

220  CALL  circle  I  XCIN)  .  0.0  .  R2IN)  »  180.0  » 

IF  I  RlIN)  -  .000001  )  240  .  240  .  230 

230  CALL  CIRCLE  I  XC I N  )  .  0.0  .  R|IN)  »  180.0  - 

240  GO  To  300 

250  CALL  PLOT  I  XCIN)  .  0.0  .1.2  ) 

CALL  PLOT  I  XCIN)  .  R2 I N )  .1.1) 

CALL  PLOT  I  Rl IN)  .  R2IN)  .  I  .  I  ) 

CALL  PLOT  I  RIIN)  .  0.0  .1.1) 

GO  TO  300 

260  IF  I  R|IN)  -  .000001  )  280  .  280  .  270 

270  call  circle  I  XCIN)  .  0.0  »  RIIN)  ,  0.0  . 

28Q  CALL  CIRCLE  I  XCIN)  .  0.0  »  R2IN)  .  PHIIN) 

300  CONTINUE 

CALL  PLOT  I  0.0  .  0.0  .1.1  ) 

CALL  FINISH  I  30  .  TITLE  ) 

END  FILE  8 

return 

END 


.  90 


-  PHI  IN)  .  -I  ) 
PHI  IN)  »  PHI  IN)  . 


PHIIN)  .  -I  ) 
.-PHIIN)  ,  I  ) 


♦LABEL 
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CDELTHA  COMPUTE  DELTA  THETAS  FOR  GEODESIC  I 

SUBROUTINE  DELTHA  (  I  ) 

dimension  R(  I  000)  .Z(  1000) ,aK(  I  000) .F(  I  000) .XI  ( I  000) .TlTLE(  12 ) . 

I RO(  I  00)  .ZO(  I  00) .ALPhaOI I  00) .C0NS( I  00 ) . W ( I  00 )  . D (  I  DO )  .ThI CK (  100)  . 
2RT( 100) .ZT( 100) .Nc( 100) .DTHETaI IDOD) .FLNgTHI 1000) 

COMMON  M.R.Z.AK.F.XI .NOGEOD .RO . ZO , ALPHAO, CONS .W » D . THI CK .NC . RT ,ZT , 
ISMaX.RMaX.ZMAX.ThMAX.JJ.T iTLE.Pi .DTHETA.fLNGTH.TSUM.FLSUM.NhIGH. 

2  NLOW  »  DISTRT  ,  ADVNCE  .  SHAFT  I  ,  SHAFT2 
CONV  t(  180.0  /  PI 

CGNS(I)  tt  RO(I)  ■>:-  SInF  (  ALPHAO(I)  /  CONV  ) 

DO  20  N  2  ,  M 


IF  1  ZOI I )  -  ZIN)  ) 

40  . 

30  . 

20 

20 

CONTINUE 

30 

I F  1  ALPHAO II)  -  90.0 

)  38 

.  33 

.  3! 

33 

IF  1  AKIN)  )  40  .  40 

.  38 

38 

N  »  N  +  1 

40 

NSTaRT  #  n  _  I 

IF  1  CONSII )  -  R  1  1  )  ) 

1  50  . 

45 

.  4! 

45 

J  tt  NSTART  +  1 

50 

J  #  J  -  1 

IF  1  Rl J)  -  CONSI  I  )  ) 

60  » 

60  . 

50 

60 

NLOW  tt  J 

IF  1  CONSI I )  -  RIM)  ) 

150  . 

70  . 

70 

70 

J  ft  NSTART 

80 

J  #  J  +  1 

IF  1  Rl J)  -  CONSI I )  ) 

90 

»  90 

.  80 

90 

NHIGH  #  J  -  1 

FOK  ft  F(NLOW)  /  AK(NLOW) 

NL  I  ft  NLOW  +  I 

term  ft  (  R(NL|)  /  CONS(l)  )**2  -  I.D 

IF  (  Term  )  92  .  92  .  95 

92  ASEC2  ft  0.0 
GO  TO  98 

95  ASEC2  ft  ATANF  (  SORTf  (  TERM  )  ) 

98  DBETA  »  ASEC2 

DTHETAI  NLOW)  ft  FOK  ♦  DBETA  *  CONV 
FLNGTH(NLOW)  ft  R(NL|)  *  FOK  *  SINF(  DBETA  ) 

NHt  tt  NHIGH  -  I 

IF  (  NH I  -  NL I  )  1 35  .  I  00  .  I  DO 

I  DO  DO  130  N  #  NL|  ,  NHI 

IF  (  AKIN)  )  120  .  f 10  .  120 

I  ID  DBETA#C0NS(  I  )*(Z(N+I  )-Z(N)  )  /  (  R  (  N  )  *SQRTF  (  R  I  N  )  «-*2  -  C0NS(I)#*2)) 
DTHETA(N)  »  DBETA  *  CONV 

FLNGTH(N)  ft  SQRTFII  Z  (  N+ I  )-Z  (  N )  )  **2  +  I  R  I  N  )  *DBETA  )  **2  ) 

GO  TO  130 

I2D  FOK  ft  ABSFI  FIN)  /  AKIN)  ) 

ASECI  ft  ASEC2 

TERM  tt  I  RIN+I)  /  CONSII)  )**2  -  I.D 

IF  I  Term  )  122  .  122  .  125 

I  22  ASEC2  #  D.O 
GO  TO  126 

125  ASEC2  tt  ATANF  I  SQRTF  I  TERM  )  ) 

126  DBETa  ft  ABSF  (  ASEC2  -  ASecI  ) 

DTHETAIN)  ft  FOK  *  DBETA  *  CONV 
RN2  tt  RIN)  »  FOK 


b 


n  n  o  n  o 


RN3  ff  R(N+|  )  *  FOK 

Term  ff  RN2*«'2  +  RN3*»2  -  2.0  *  RN2  *  RN3  *  COSF  (  DBETA  ) 

IF  (  Term  >  127  »  i27  .  129 

I  27  FLNGTHIN)  #  n.H 
GO  TO  130 

129  FLNGTHIN)  #  SQRTF  I  TERM  ) 

I  30  CONTINUE 
135  CONTINUE 

FOK  /i  ABSF  I  F(NHIGH)  /  AK(NHIGH)  ) 

DBETA  »  ABSF  (  ASEC2  ) 

DTHETA(NHIGH)  »  FOK  *  DBETA  *  CONV 
FLNGTH(NHIGH)  »  FOK  *  RINHIGH)  *  SINF  (DBETA) 

TsuM  #  n.n 
FLSUM  #  n.n 

DO  140  N  ^  NLOW  *  NHIGH 
TSUM  »  TSUM  +  DTHETA(N) 

140  FLSUM  tf  FLSUM  +  FLNGTH(N) 

TSUM  ff  2.0  *  TSUM 
FLSUM  tf  2.0  *  FLSUM 
GO  TO  160 

150  CONS(I)  tt  MAXIF  (  R(|)  .  RIM)  ) 

ALPHAOII)  #  ATaNFI  CONS! I)/  SQRTFI  RO(I)**2  -  CONS(I)**2)  )  *  CONV 
WRITE  OUTPUT  TaPE  6  ,  8000  .  I  .  ALPHAOII) 

8000  FORMAT  I  33H0  TURN-AROUND  RADIUS  FOR  GEODESIC  .13  .  56H  IS  LESS  T 
IHAN  R(|)  OR  RIM)  -  STARTING  ANGLE  CHANGED  TO  .  FID. 6  ) 

GO  TO  45 
160  CONTINUE 
RETURN 
END 

»LABEL 


CNOCIRC  COMPUTE  NUMBER  OF  CIRCUITS  TO  GIVE  THICKNESS 

SUBROUTINE  NOCIRC  I  I  .  KK  »  LA  ) 

I  IS  GEODESIC  NUMBER 

KK  IS  option  IN  determining  NUMBER  OF  CIRCUITS  PER  PATTERN 
LA  IS  option  TO  ADJUST  STARTING  ANGLE  OR  DISTORT  GEODESIC 

DIMENSION  R(  I  000  )  .Z-l  I  000  )  .  aKI  I  000)  »F  U  000  )  .X  I  I  I  000  )  .TITLE  I  12). 

I ROI  I  DO)  .ZOI  I  00) .ALPHaOI  I  00) .CONS  I  I  00 ) . W ( I  00 )  .  D I  I  00  )  .  Th  I CK I  I  DO  )  . 
2RT(  IDO)  .ZTI  100)  .NC(  100) .DTHETaI  1000) .FLNgTHI  I  ODD) 

COMMON  M.R.Z.AK.F.XI  .NOgEOD »R0 ,Z0 , ALPhAO.CONS  .W . D . TH I CK  ,NC . RT .ZT. 
ISMAX.RMAX.ZMAX.THMAX.JJ.T iTLE ,P I . DTHETA  .FLNGTH . TSUM . FLSUM, NHI GH . 

2  NLOW  ,  DISTRT  ,  ADVNCE  ,  SHAFT  I  ,  SHAFT2 
PGLASS  #  I  .0 

IF  I  D( I )  )  10.10.20 

10  Dll)  tf  .00  1 

write  output  tape  6,1400  .  I 

1400  FORMAT  I  35H0  DIAMETER  OF  ROVING  FOR  GEODESIC  ,  I3  .  23H  NOT  GIVE 
IN  -  .001  USED  ) 

20  IF  I  W( I )  )  30  ,  30  .  40 

30  W( I )  #  .1 

WRITE  OUTPUT  TAPE  6  ,1410  .  I 

1410  FORMAT  I  32HD  WIDTH  OF  ROVING  FOR  GEODESIC  ,  13  .  2IH  NOT  GIVEN 
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I-  .1  USED  ) 

40  IF  (  THICK(I)  )  42  ,  42  ,  48 

42  THICK(  I  )  tt  2.0  *  DU) 

WRITE  OUTPUT  TAPE  6.1415.  I 

1415  FORMAT(32H0  DESIRED  THICKNESS  FOR  GEODESIC  ,  I3.  29H  NOT  SPECIFIED 
I  .  2  D( I )  USED  ) 

48  IF  (  ZO(I)  )  50  .  50  .  120 

50  IF  (  RO(I)  )  60  .  60  .  70 

60  write  output  Tape  6. 142a  .  i 

1420  FORMAT  (  46HD  STARTING  STATION  NOT  SPECIFIED  FOR  GEODESIC  ,  13  , 

I  I6H  CANNOT  COMPUTE  ) 

GO  TO  500 

70  DO  80  N  #  2  .  M 

IF  (  RO( I )  -  R(N)  )  90  .  85  .  80 

80  CONTINUE 

WRITE  OUTPUT  TAPE  6.1430  .  I 

1430  FORMAT  (  72HD  WITH  ZO  NOT  GIVEN  .  COULD  NOT  DETERMINE  STARTING  SE 
ICTION  for  GEODESIC  .  13  .  I  OH  USING  RO  ) 

GO  TO  500 
85  N  »  N+l 
90  NSTART  #  N  -  I 

IF  <  AK(NSTaRT)  )  100  .  MO.  100 

too  ZO(I)  »  (  RO(I)  -  R(NSTaRT)  )  /  AK(NSTART)  +  Z(NSTART) 

GO  TO  160 

MO  ZO(  I  )  H  Z(NSTART) 

GO  TO  160 

120  DO  130  N  »  2  .  M 

IF  (  ZO(I)  -  Z(N)  )  140.  135.  130 

I  30  CONTINUE 

WRITE  OUTPUT  TAPE  6.1440.1 

1440  FORMAT  (  63H0  USING  Zo  .  COULD  NOT  DETERMINE  STARTING  SECTION  FOR 
I  GEODESIC  .  13  ) 

GO  To  500 
135  N  N  +  I 
140  NSTART  #  N  -  I 

IF  (  ROM)  )  150  .  150  »  160 

150  RO(I)  »  AK(NSTaRT)  *  (  ZO(I)  -  Z(NSTaRT)  )  +  R(NSTART) 

160  CONTINUE 

IF(  ZT( I)  )  170  .  170  .  240 
170  IF(  RT(I)  )  180  .  180  .  190 

180  write  output  Tape  6  .1450  .  i 

1450  FORMaT(3IH0(RT,ZT)  NOT  GIVEN  FOR  GEODESIC  .  13  . I 3H  (RO.ZO)  USED  ) 
185  RT(  I  )  ft  RO(  I  ) 

ZT(  I  )  ZO(  I  ) 

NTH  »  NSTART 
GO  TO  280 

190  DO  200  N  #  2  .  M 

IF  (  RT(I)  -  R(N)  )  210.  205.  200 

200  CONTINUE 

write  output  Tape  6.i46o  .  i 

1460  format  (  83H0  WITH  ZT  NOT  GIVEN  .  COULD  NOT  DETERMINE  SECTION  TO  C 
lOMPUTE  THICKNESS  FOR  GEODESIC  .  13  ,  39H  USING  RT  ,  SO  (RO.ZO)  US 

2Eti  FOR  (RT.ZT)  ) 

GO  TO  185 
205  N  »  N+l 
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2  I  0  NTH  #  N  -  I 

IF  (  AK(NTH)  )  220  .  230  .  220 

220  ZT(I)  #  (  RT(I)  -  R(NTH)  )  /  AK(NTH)  +  Z(NTH) 

GO  To  280 
230  ZT(  I  )  ff  Z(NTH) 

GO  TO  280 

2A0  DO  250  N  i¥  2,  M 

IF  (  ZT(I)  -  Z(N)  )  260  »  255,  250 

250  CONTINUE 

WRITF  OUTPUT  TAPE  6,1470  ,  I 

1470  format  (  80H0  USING  GIVEN  ZT  ,  COULD  NOT  DETERMINE  SECTION  TO  COMP 
|UTE  THICKNESS  FOR  GEODES  I C , I  3  ,28H  So  (RO,ZO)  USED  FOR  (RT,ZT)  ) 

GO  TO  185 
255  N  #  N  +  I 
260  NTH  »  N  -  I 

IF  (  RT(I)  )  270  ,  270  ,  280 

270  RT(I)  #  AK(NTH)  *  (  ZT(I)  -  Z(NTH)  )  +  R(NTH) 

280  CONTINUE 

290  CONS(I)  #  RO(n  *  SINF  (  ALPHAOI  1)  *  PI  /  180.0  ) 
call  PeRcOV  (  RT(I)  ,  nth  ,  I  ,  PERCNT  ,  HANGL  ) 

IF  (  PERCNT  )  300  ,  300  ,  310 

300  RT(  I)  If  CONS(  I  ) 

ZT(  I  )  ft  0.0 

WRITE  OUTPUT  TAPE  6  ,1480,  I 

1480  format  (  42H0  COVERAGE  AT  (RT,ZT)  IS  ZERO  FOR  GEODESIC  ,  13  ,  38H 
|,  TURNAROUND  POINT  USED  FOR  (RT,ZT)  ) 

GO  TO  160 

310  IF  (  KK  )  320  ,  320  ,  330 

320  B  ft  THICK(I)  /  (  Dd)  *  PERCNT  )  *  PGLASS 
GO  TO  340 

330  B  ft  2.0  /  PERCNT 
340  NB  #  B  +  .5 

CALL  DELTHA  (  I  ) 

NLOW  ft  NLOW 
NHIGH  ft  NHIGH 
RVN  ft  TSUM  /  360.0 
INTGR  ft  RVn 

FRACT  ft  RVN  -  FLOATF(  INTGR  ) 

A  #  B  *  FRACT 
NA  #  A  +  .5 

IF  (  NA  )  350  ,  350  ,  360 

350  NA  ft  I 

GO  TO  388 

360  IF  (  NA  -  NB  )  380  ,  370  ,  375 

370  NA  #  NB  -  I 
GO  TO  388 

375  INTGR  ft  INTGR  +  I 
NA  #  NA  -  NB 

380  CALL  NOFACT  (  NA  ,  NB  ) 

388  CONTINUE 

EPS  ft  .000001 

IF  (  LA  )  382  ,  382  ,  390 

382  CONTINUE 

NAA  #  NB  *  INTGR  +  NA 

call  adjust  (  I  ,  NAA  ,  NB  ,  RVN  ,  EPS  ,  LL  ) 


IF  (  LL  )  384  ,  384  ,  39n 

384  CALL  PERCOV  (  RT(I)  .  NTH  .  I  .  PERCNT  ,  HANGL  ) 

DISTRT  U  I  .0 
GO  TO  420 

390  RVN2  n  FLOaTFI  INTGR  )  +  FLOATFI  NA  )  /  FLOATF  (  NB  ) 

DISTRT  »  RVN2  /  RVN 
DO  410  N  #  NLOW  *  NHI6H 
410  DTHETA(N)  H  DTHETA(N)  *  DISTRT 
TSUM  #  TSUM  *  DISTRT 

420  AN  H  ThICK(I)  /  (  Dd)  *  PERCNT  )  *  PGLASS 
NOPATN  #  AN  /  FLOATFI  NB  )  +  .5 

NCPERP  #  NB 
NCI  I )  n  NCPERP  *  NOPATn 

THNESS  n  FLOATF  I  NC 1 1 )  )  *  Dll)  *  PERCNT  /  PGLASS 
ZLOW  #  I  cONSIi)  -  RINLOW)  )  /  AKINLOW)  +  ZINLOW) 

ZHIGH  #1  CONSII)  -  RINHIGH))  /  AKINHIGH)  +  ZINHIGH) 

write  output  Tape  6  ,1490  ♦  i»  ncid  .nopatn  ,  ncperp  .  intgr» 

I  NA»  NB»  ThNESS  ,  CONSIi)  ,  ZLOW  »  ZHIGH 
1490  format  I  lOOHQ  NO.  OF  NO.  OF  CIRC.  PER  RATIO  OF 

I  INTEGERS  thickness  TURNAROUND  STATIONS  /  I D6H  GEODESIC  Cl 

2RCUITS  patterns  PATTERN  N  +  A  /  B  AT  IRT.ZT)  RA 

3DIUS  Z  lower  Z  upper  /  IHO  .  14.  6110  .  4FID.6  ) 

write  output  Tape  s.isoo.  distrt 
1500  format  I  23H0  DISTORTION  FACTOR  »  ,  FlO.6  ) 

write  output  Tape  6. i5io»  in.dthetain)  .  flngthin).  m  nlow.nhigh  ) 

1510  FORMAT  I  39HaSECTlON  DELTA  THETA  FILAMENT  LENGTH  /  I  IH  »  14  , 

I  2FI6.6  )  ) 

WRITE  OUTPUT  TAPE  6  ,1520  ,  TSUM  ,  FLSUM 
1520  FORMAT  I  8H0CIRCUIT  ,  FI3.6  .  FI6.6  /  IHO  ) 

500  CONTINUE 
RETURN 
END 

*LABEL 


CGEOPLT  plot  geodesic  ON  DEVELOPED  SURFACE 

SUBROUT  INE  GEOPLT  I  I  ,  NUM  ) 

dimension  RI 1000) .ZI 1000) .A*^( InOO) .F I  1000) .XI  I  1000) .TiTLeI  12) . 

I ROI  I  00 )  .zoi  I  00 ) .ALPHAOI  I  00) .cons I  I  00 ) .W  I  I  00)  .Dl  I  00)  .ThICKI  I  00  )  , 
2RT  I  I  00)  .ZTI  I  00 ) ,NC( I  00) .DTHETaI  I  000 ) .FLNgTHI  1000) 

DIMENSION  RlllOO  )  .  R2 I  I  00  )  .  PHUlOD  )  .  XCIlOO  ) 

COMMON  M.R.Z.AK.F.XI  .NOGEOD .RO ,ZO , ALPHAO , C ONS , W , D  .THI CK  ,NC . RT ,Z T , 
ISMAX.RMAX.ZMAX.THMAX.JJ.TITLE.PI .DTHETA.FLNGTH.TSUM.FLSUM.NHIGH. 

2  NLOW  ,  DISTRT  ,  ADVNCE  .  SHAFT  I  ,  SHAFT2 
COMMON  AA ,BB.CC,DEL,DELRH0,NSTART 
COMMON  LLL  ,  RHOMIN  ,  FR  ,  TMIN 
COMMON  R|  ,  R2  ,  PHI  ,  XC 
KOUNT  «  0 

CONV  #  PI  /  180.0 

THETa  #  0.0 

write  output  Tape  0.4010  .  i 

4010  format  I  28H  DEVELOPED  PLOT  OF  GEODESIC  ,  13  .  5H$  ) 

READ  INPUT  TAPE  0,4020.  I  TITLEIK)  ,  K  »  7,12  ) 

4020  FORMAT  I  6A6  ) 

CALL  setup  I  TITLE  ) 


QOOOGI  )  20  *  20  . 


DO  10  N  »  I  >  M 
IF  (  RMaX  -  R(N)  -  . 


1  0 

CONTINUE 

20 

NSTART  ft  N 

ZZ  ft  Z(  NSTART) 
UPDOWN  #  1 .0 

30 

IE  (  AKIN)  ) 

40 

♦  42  0» 

70 

40 

AKK  »  -  1 .0 

IE  (  UPDOWN  ) 

50 

*  50  . 

60 

50 

RR  ti  R|  IN) 

GO  TO  too 

60 

RR  ft  R2  (N  ) 

GO  TO  100 

70 

AKK  #  1.0 

IE  (  UPDOWN  ) 

80 

»  80  » 

90 

80 

RR  tt  R2  IN) 

GO  TO  100 

90 

RR  ft  R  1  IN  ) 

1  00 

AOF  tt  ABSf  I  AKIN) 

/  FIN) 

) 

BETA  #  AOF  *  theta 

XO  #  AKK  *  RR  *  COSE  (  BETA  *  CONV  ) 
YO  »  RR  *  SINE  (  BETA  *  CONV  ) 

IE  (  AKK  )  I  10  »  I  10  »  160 


1  1  0 

IF 

1  N  -  NHIGH 

) 

130  ♦ 

170  »  170 

1  30 

IF 

1  UPDOWN  ) 

140 

140  » 

1  50 

1  40 

RE 

tt  R2  IN) 

GO 

TO  2  10 

1  50 

RE 

tt  R  1  IN  ) 

GO 

TO  210 

1  60 

IF 

1  N  -  NLOW  ) 

170 

•  170 

»  180 

1  70 

RE 

tt  RR 

BETAZ  tt  beta  + 

2.0 

* 

AOF  * 

DTHETAIN) 

GO 

TO  220 

180 

IE 

1  UPDOWN  ) 

190 

♦ 

190  » 

200 

1  90 

RE 

tt  R  1  IN) 

GO 

TO  2  10 

200 

RE 

tt  R2  IN) 

210 

BETAZ  tt  BETA  + 

AOF 

* 

DTHETAIN) 

220 

CONT INUE 

XD  »  AKK  *  RE  *  COSE  (  BETAZ  *  cONV  ) 
YD  RE  *  SINE  (  RETaZ  *  CONV  ) 


IE  (  BETAZ  -  PHI(N)  )  230  »  360  »  360 

230  XORE  #  XO  +  XC(N) 

XDRE  #  XD  +  XC(N) 

CALL  PLOT  (  XORE  »  YO  ,  I  ,  2  ) 

CALL  PLOT  (  XDRE  .  YD  .  I  *  I  ) 

theta  ft  BETAZ  /  AOE 
IE  (  AKK  )  240  »  290  »  290 

240  IE  (  N  -  NHIGH  )  250  .  270  *  270 
250  IE  (  UPDOWN  )  280  »  280»  260 
260  ZZ  »  Z(N+I ) 

N  #  N  +  I 
GO  TO  340 
270  UPDOWN  #  -  I .0 
280  ZZ  ft  Z(  N) 
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N  #  N  -  I 
GO  TO  340 

290  IF  (  N  -  NLOW  )  320  ♦  320  »  300 
300  IF  (  UPDOWM  )  310  »  310  ♦  330 
310  ZZ  »  Z ( N) 

N  »  N  -  I 
GO  TO  340 
320  UPDOWN  n  I .0 
330  ZZ  »  Z(N+I ) 

N  #  N  +  I 
340  continue 

IF  (  ZZ  -  Z(NSTaRT)  )  30  »  350  ♦  30 

350  KOUNT  #  KOUNT  +  | 

IF  (  KOUNT  -  2  *  NUM  )  30  »  510  .  510 

360  Al  U  SiNF  (  PHKN)  *  CONV  ) 

Bl  »  -  AKK  *  COSF  (  PHI(N)  *  CONV  ) 

IF  (  ABSf  (  XO  -  XD  )  -  .0001  )  370  »  370  *  380 

370  A2  #  1.0 
B2  #  0.0 
D2  #  XO 
GO  TO  390 

380  SLPE  #  (  YD  -  YO  )  /  (  XD  -  XO  ) 

A2  »  -  SLPE 
B2  »  1.0 

D2  »  YO  -  SLPE  ♦  XO 
390  DENOM  #  A I  »  B2  -  A2  *  B I 

IF  (  ABSF  (  DENOM  )  -  .0001  )  410  »  410  »  400 

400  XI  #  (  -  Bl  »  D2  )  /  DENOM 

YI  «  Al  *  D2  /  DENOM 
XORE  »  XO  +  XC(N) 

XIRE  «  XI  +  XC(N) 

CALL  PLOT  (  XORE  .  YO  .  I  .  2  ) 

CALL  PLOT  (  XIRE  .  YI  ,  I  ,  I  ) 

XO  «  AKK  »  SQRTE  (  Xl»*2  +  YI*«2  ) 

YO  #  0.0 

BETAZ  »  BETAZ  -  PHI (N) 

GO  To  220 

410  write  output  Tape  6  ,  40aO  *  N.aKCN) ,AI  »BI  »A2»b2»D2»PHI (N) ,X0,Y0, 

I  XD,  YD 

4000  EORMAT  (  60HI  LINE  CONNECTING  (XO.YO)  AND  (XD,YD)  IS  PARALLEL  TO  L 
I INE  2  /  IHO  ,  13  *  I IFIO.4  ) 

GO  TO  510 

420  IF  (  UPDOWN  )  430  ,  430  ,  440 

430  XO  #  R I (N ) 

XD  #  XC(N) 

ZZ  #  Z(N) 

NN  #  N  -  I 
Go  TO  450 
440  XO  #  XC(N) 

XD  #  Rl (N) 

ZZ  »  Z( N+l ) 

NN  #  N  +  I 

450  YO  #  Theta  *  R(N)  *  conv 

YD  YO  +  R(N)  *  DTHETA(N)  *  CONV 
SLPE  »  (  YD  -  YO  )  /  (  XD  -  XO  ) 
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460  continue 

IF  (  YD  -  R2(N)  )  4RG  »  480  »  470 

470  YI  ff  R2 (N) 

XI  S'  (  SURE  *  XO  +  YI  -  YO  )  /  SLPE 

CALL  PLOT  (  XO  ,  YO  ,  I  ,  2  ) 

CALL  PLOT  (  XI  ,  YI  ,  I  ,  I  ) 

XO  ^  XI 

YO  »  0.0 

yd  #  Yd  -  R2(N) 

GO  TO  460 

480  CALL  PLOT  (  XO  ,  YO  >  I  ,  2  ) 

CALL  PLOT  (  XD  ,  YD  ,  I  ,  I  ) 

THETa  i¥  Yd  /  (  R(M)  *  CONV  ) 

IF  (  ZZ  -  Z(NSTART)  )  500  .  490  .  500 

490  KOUNT  ff  KOUnT  +  | 

IF  (  KOUNT  -  2  *  NUM  )  500  »  510  »  510 

500  N  »  NN 

GO  TO  30 
510  CONTINUE 

CALL  FINISH  (  30, TITLE  ) 

END  FILE  8 

RETURN 

FND 

»LABFL 


CDIVSUR  DIVIDE  UP  SURFACE 

subroutine  DIVSUR  (  STEP  ) 

DIMENSION  R( I  000 ) ,Z( I  000 ) ,aK(  I  000) »F( 1000 ) ,XI  ( I  000 ) » T i T LE t  12), 

I RO ( I  00 ) »ZO ( I  00 ) , ALPHaOC I  00) ,CONS( I  00) ,W( I  00)  ,D(  I  00)  ,ThI CK(  I  DO  )  , 
2RT ( I  00)  ,ZT( I  00 ) ,NC(  I  00 ) ,DTHETa( I  000) ,FLNgTh(  1000) 

DIMENSION  MANGLE (  I  00) ,ThNESS( | 00 ) 

COMMON  M,R,Z,AK,f,XI ,NOgEOD ,R0 ,Z0 , aLPHAO, CONS , W , D , THI CK  ,Nc , RT , ZT , 
|SMAX,RMAX,ZMAX,THMAX,JJ,TlTLE,Pl jDTHETA  ,FLNGTH , TSUM  ,  FLSUM , Nh I GH , 

2  NLOW  ,  DISTRT  ,  ADVNCE  ,  SHAFT  I  ,  SHAFT2 
ThMaX  #  0.0 
J  #  0 
S  #  0.0 
MM  I  #  M  -  I 
DO  100  N  I  ,  MMI 
IF  (  AKIN)  )  5  ,  50  ,  5 

5  S  #  X  I  ( N) 

AKOFN  #  AKIN)  /  FIN) 

RPX  »  RIN)  -  AKOFN  *  XI  IN) 

ZPX  #  ZIN)  -  XI  IN)  /  FIN) 

I  0  J  #  J  +  I 

RR  ff  RPX  +  AKOFN  *  S 
ZZ  #  ZPX  +  S  /  FIN) 

15  SUMTH  #  0.0 

DO  20  I  »  I  ,  NOGEOD 

CALL  PERCOV  I  RR  ,  N  ,  I  ,  PERCNT  ,  HANGL  ) 

MANGLE!  I  )  ff  HANGL 

THNESSII)  ff  FLOATF  I  NCII)  )  *  DID  *  PERCNT 
SUMTh  a-  SUMTH  +  THNESSli) 

20  CONTINUE 


r 


THMAX  »  MAX  IF  (  THMAX  .  SUMTH  ) 

RFINAL  #  RR  +  <SUMTH  /  F(N) 

ZFINAL  ft  11  -  SUMTH  *  AKOFN 

WRITE  TAPE  I  »  S  ,  RR  »  ZZ  »  (  HANGLE(I)  »  I  I  »  NOGEOD  )  » 
I  (THNESS(I)  ,  I  I  ,  NOGEOD  )  ,  SUMTH  ,  RfiNAL  ,  ZFINAL 
IF  (  S  -  XI(N+I)  +  STEP  )  25  ,  30  »  3D 
25  S  #  S  +  STEP 
GO  TO  ID 

3D  IF  (  S  -  XKN+I)  +  .000001  )  35  ,  IDO  .  IDO 

3.5  S  ft  X  I  (  N+  I  ) 

RR  ft  R(N+|  ) 

11  ft  Z(N+I  ) 

J  #  J  +  I 
GO  TO  15 
50  J  »  J  +  I 
S  ft  X  I  (N) 

RR  #  R ( N) 

ZZ  ft  Z  (  N  ) 

SUMTh  ft  0.0 

DO  80  I  I  ,  NOGEOD 

CALL  PERCOV  (  RR  .  N  ,  I  ,  PERCNT  ,  HANGL  ) 

HANGLEI  I  )  ft  HANGL 

THNESS(I)  #  D(I)  »  PERCNT  *  FLOATF  (  NC(I)  ) 

SUMTH  ft  SUMTH  +  THNESS(I) 

80  CONTINUE 

THMAX  ft  MAX  IF  (  THMAX  .  SUMTH  ) 

RFINAL  #  RR  +  SUMTH 
ZFINAL  ft  11 
J2  H  I 

90  write  Tape  i  »  s  ,  rr  ,  zz  ♦  (  hangle(I)  »  i  »  i  »  nogeod  i  » 

I  (  THNESS(I),  I  #  I  ,  NOGEOD  )  ,  SUMTH  .  RFINAL  .  ZFINAL 
IF  (  J2  -  2  )  95  .  ino  ♦  ino 

95  S  ft  X  I  (  N+  I  ) 

J  #  J  +  I 
ZZ  ft  Z  (  N+  I  ) 

ZFINAL  ft  11 
J2  #  2 
GO  TO  90 
100  CONTINUE 
JJ  #  J 
SMAX  ft  S 
END  FILE  I 

return 

END 

♦label 


crzplot  plot  r  and  z  versus  s 

SUBROUTINE  RZPLOT  (  XO  ,  YO  ,  XL  ,  YL  ) 

dimension  R(  I, ODD)  .Z(  1000)  »aK(  I  000)  .Ft  I  ODD)  »XI  (  I  ODD)  .TiTLeI  12). 
|R0( 100) .Z0( 100) .ALPhaO( I  DO) .CONSI I  DO ) . W ( I EQ ) . D(  I  00 ) .Th  I CK  (  I  DO )  . 
2RT ( IDO) .ZT ( 100) .NC( 100) .DThETa ( I  000 ) >FLNgTh(  I  ODD ) 

DIMENSION  XAX(I2)  ,  A(I2)  ,  Y(2I0) 

COMMON  M.R.Z.AK.F.XI  .NOGEOD. RO.ZO .ALPHAO , CONS .W . D » TH I CK .NC » RT ,ZT » 
I SMAX. RMAX, ZMAX. THMAX. JJ. TITLE. PI .DTHETA .FLNGTH . TSUM . FLSUM . NH I GH . 


COCQCDCDCDCDCOCD 


2  NLOW  ,  DISTRT  ,  ADVNCE  »  SHAFT  I  »  SHAFT2 
TITLE(7)  ff  606060605161 
TITLE(8)  »  514421676060 
TITLE(9 )  #  2  I  452460716  1 
TITLE(  IO)iV  714421676060 
TITLE(ll)i!(  536060606060 
XAX(I)  #  606060606261 

XAX(2)  «  624421676060 

XAX(3)  #  536060606060 

DIVX  )f  10.0  /  XL 

DIVY  #  10.0  /  YL 

CALL  Setup  (  title  ) 

call  AXPLOT  (  XO.yO.xL  ,  yl  .DIVX.DIVY,  .1,.I.  5,5,50»XAX  » 

I  TITLE(7)  ) 

INK  »  2 

DO  3  N  #  I  »  M 
Yl  #  R(N)  *  YL  /  RMAX  +  YO 
SI  #  X  I  (N )  *  XL  /  SMAX  +  XO 
CALL  PLOT  (  SI  .  Yl  ,  I  ,  INK  ) 

INK  #  I 

3  CONTINUE 
INK  »  2 

DO  4  N  »  I  »  M 
Yl  »  Z(N)  *  YL  /  ZMAX  +  YO 
SI  #  X I  (N )  *  XL  /  SMAX  +  XO 
CALL  PLOT  (  SI  .  Yl  .  I , INK  ) 

INK  #  I 

4  CONTINUE 

write  output  Tape  o»40io»  xax(2)  *  smax 
4010  format  (  A6  .  IH#  .  F8.4  .  I H$  ) 

READ  INPUT  TAPE  0.4020.  (A(K).  K  #  1.3  ) 

4020  FORMAT  (  3A6  ) 

XX  #  XL . /  4.0  +  XO 
YY  »  YL  +  .5  +  YO 

CALL  LETTER  (  XX  .  YY  .  50  .  52  .  A  ) 

write  output  Tape  0.4010.  title(8).  rmax 

READ  INPUT  TAPE  0.4020.  (  A(K)  ,  K»  I ,  3  ) 

YY  #  YL  +  .3  +  YO 

CALL  letter  (  XX.  YY ,  50  .  52  .  A  ) 

write  output  Tape  o.aoid  .  titleUO)  .  zmax 
read  input  tape  0.4020. (  A(K)  ,  K  #  I  ,  3  ) 

YY  #  YL  +  . I  +  YO 

CALL  LETTER  (  XX.  YY ,  50.  52  .  A  ) 

YY  #  YO  +  YL  +  2.0 

call  letter  (  XO.  YY,  50.  52  .  TITLE  ) 

CALL  FINISH  (  30.  TITLE  ) 

END  FILE  8 

return 

END 

♦LABEL 


canglpl  plot  helix  angle  versus  S 

subroutine  ANGLPL  (  XO  ,  YO  ,  XL  ,  yl  ) 

dimension  R(  I  000)  .Z(  1000)  .AX(  I  000) .F( I  000)  .X  I  (  I  000) .TiTLe(  12). 


□D  CD  CD  CD  CD  CD 


I RO (  I  00 ) »Z0( I  00 ) .ALPHA0( I  00) »CONS( I00)»W(|00)>D(I00)»THICK(I00)» 
2RT( 100) >ZT( 100) .Nc( 100) .DTHETa( 1000) »FLNgTH( 1000) 

DIMENSION  XAX(I2)  »  A(I2)  ,  Y(2lD) 

COMMON  M.R,Z,AK,F.XI ,N0GE0D»R0 ,ZO , ALPHAO ,CONS . W » D . TH I CK . NC . RT ,ZT 
|SmAX,RmAX,ZmAX,THmAX,JJ»TITlE,PI »DTHETA,FLNGTH,TSUm»FL5UM.NHIGH, 
2  NLOW  ,  DISTRT  ,  ADVnCE  >  SHAFT  I  >  SHAFT2 
TITLE!?)  n  603025433167 
TITLE(8)  »  602145274325 

TITLE(9)  #  536060606060 

XAX(I)  U  606060606261 
XAX(2)  n  624421676060 
XAX(3)  n  536060606060 
DIVX  #  10.0  /  XL 
DIVY  n  9.0  /  YL 
CALL  SETUP  (  TITLE  ) 

CALL  AXPLOT  (  XO . YO.XL , YL .D I VX » DI VY , . I  ,10.0  ,  5,  3  ,  50  ,  XaX  , 
I  TITLE!?)  ) 

DO  20  I  #  I  ,  N06E0D 
INK  #  2 
REWIND  I 

DO  20  J  »  I  ,  JJ 
12  »  I  +  2 

READ  TAPE  I  ,  S  ,  !  Y!K)  ,  K  #  I  ,  12  ) 

IF  !  Y! 12)  )  20  ,  20  ,5 
5  YI  #  Y! 12)  *  YL  /  90.0  +  YO 
SI  #  S  *  XL  /  SMAX  +  XO 
CALL  PLOT  !  SI  ,  Y1  ,  I  ,  INK  ) 

INK  It  I 
20  CONTINUE 

write  output  Tape  o  ,  4oio  ,  xax!2)  ,  smax 

4010  FORMAT  !  A6  ,  IH#  ,  F8.4  ,  IH$  ) 

READ  INPUT  TAPE  0,4020,  !  AlK),  K  »  1,3  ) 

4020  FORMAT  I  3A6  ) 

XX  XO  +  XL  /  4.0 
YY  #  YO  +  YL  +  .  I 

CALL  letter  !  XX  ,  YY  ,  50  ,  52  ,  a  ) 

YY  #  YO  +  YL  +  2.0 

CALL  Letter  !  xo  ,  yy  ,  5n,  52  ,  title  ) 

CALL  FINISH  !  30,  TITLE  ) 

END  FILE  8 
REWIND  I 
RETURN 
END 

♦LABEL 


CTHPlOT  PLOT  thickness  VS  S  FOR  GEODESIC  I 

subroutine  THPLOT  ! I  ,  xo  ,  YO  ,  XL  ,  YL  ) 

dimension  RI  I  000 ) ,Z ! 1000 ) ,AK!  I  000) ,F !  I  000 ) ,X I  !  I  000) ,TlTLE!  12), 

1  RO!  I  00)  ,Z0!  I  00) ,ALPHA0!  | 00) , CONS  I  I  00 ) , W  !  I  00 )  , D !  I  00  )  ,TH I CK !  I  00 )  , 
2RT !  I  00)  ,ZT!  I  00 ) ,NC!  100) ,DTHETA !  I  000) ,FLNGTH! 1000) 

DIMENSION  XAX!|2)  ,  A!|2)  ,  B!|2)  ,  Y!2I0) 

COMMON  M,R,Z,AK,F,X|  ,NOGEOD ,R0 , ZO , ALPHAO , CONS , W , D , TH I CK  ,NC , RT  ,ZT 

I SMAX, rmax,zmax,thmax,jj, Title ,pi ,dtheta,flngth,tsum,flsum,nhigh, 

2  NLOW  ,  DISTRT  ,  ADVNCE  ,  SHAFT  I  ,  SHAFT2 


CQ  CO  Qj  CD  CQ  CD  CD 
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XAX(  I  )  tf  60606062606  I 
XAX(2)  #  6D62442I676D 

XAX(3)  tf  536060606060 
TITLE(7)  U  606063606160 
TITLE(8)  tt  634421676026 
TiTLE(9)  #  465160272546 
TITLE( 10)  #  242562312360 
DIVX  «  10.0  /  XL 
DIVY  #  lO.G  /  YL 
XLD  »  . I 

YLD  #  .  I 

NX  »  5 

NY  #  5 

II#  50 

write  output  Tape  o»iqi5  .  xax(2)  .  smax 

1015  FORMAT  (  A6  ,  IH#  .  F8.4  .  I H$  ) 

READ  INPUT  TAPE  0.1020.  (  AIK).  K  »  I .3  ) 

1020  FORMAT  (  3A6  ) 

write  OUTPUT  TAPE  0.1025.  THMAX 
1025  FORMAT  (  2nHMAXIMUM  THICKNESS  »  .F9.6  .  IH$  ) 

READ  INPUT  TAPE  0.1030.  (  B(K)  .  K  #  I,  5  ) 

1030  FORMAT  (  5A6  ) 

50  WRITE  OUTPUT  TAPE  0.  1005  .  I 
1005  FORMAT  (  12  .  |H$  ) 

READ  INPUT  TAPE  O.IOIO.  TITLE(II) 

1010  FORMAT  (  A6  ) 

60  REWIND  I 

CALL  SETUP!  TITLE  ) 

CALL  AXPLOT  (XO.YO.XL.YL.DIVX.DIVY.XLD.YLD.NX.NY.il .XAX.T  I  TLe  (  7 ) ) 
INK  #  2 

DO  70  J  #  I  .  JJ 
12  #  I  +  NOGEOD  +  2 

READ  TAPE  I.  S.  (  Y(K)  .  K  #  I  .  12  ) 

Si  #  S  *  XL  /  SMAX  +  XO 
YI  #  Y( 12 )  *  YL  /  THMAX  +  YO 
CALL  Plot  (  si .  yi  .  i  .  ink  ) 

70  INK  #  I 

XX  #  XO  +  XL  /  4.0 
YY  #  YO  +  YL  +  .3 

CALL  LETTER  (  XX.  YY  .  50  .  52  .  A  ) 

YY  #  YO  +  YL  +  .1 

CALL  Letter  (  xx.  yy,  50  .  52  .  b  ) 

YY  #  YO  +  YL  +  2.0 

CALL  letter  (  XO  .  yy  ,  50  .  52  .  TITLE  ) 

CALL  FINISH  (  30  .  TITLE  ) 

80  end  file  8 

REWIND  I 

RETURN 

END 

*LABEL 


CSUMPLT  plot  sum  of  THICKNESS  VERSUS  S 

subroutine  SUMPLT  (  XO  ,  yO  ,  XL  ,  YL  ) 

DIMENSION  R(IDOO).Z(|DOO).AK(| 000) .F ( I  000)  .X |  (  I  000) .TITLE!  12). 
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I R0(  I  00) »Z0(  I  DO ) .ALPHA0( I QQ) .CONS( I  00 ) *W( I  DO ) ♦D(  I  00 ) »ThICK(  IDO). 

2RT ( I  00)  »ZT (  I  00 ) >NC(  I  00) ,DTHETa(  1000 ) .FLNGTH(  1000) 
dimension  XAX(|?)  .  A(|2)  ♦  B(|2)  .  Y(||0) 

COMMON  M»R»Z,AK»F»X| »N06E0D*R0 »Z0 . ALPHAO , CONS  * W . D . TH I CK »NC » RT ,21 , 
|SMAX,RMAX,ZMAX,THMAX,JJ,TlTLE.Pl jDTHETA *FLNGTH * TSUM . FLSUM , NHI GH * 

2  NLOW  .  DISTRT  ,  ADVNCE  »  SHAFT  I  ,  SHAFT2 
B  XAX( I )  n  606060626061 

B  XAX(2)  it  606244216760 

B  XAX(3)  it  536060606060 

B  TITLE(7)  it  606063606160 

B  TiTLElS)  #  634421676026 

B  TITLE(9)  it  465160272546 

B  TITLEIIO)  it  242562312360 

B  55  TITLE!  II)  it  606264445360 
DI VX  #  I  0.0  /  XL 
DIVY  it  10.0  /  YL 
XLD  it  .  I 
YLD  it  .  I 
NX  »  5 

NY  #  5 

\l  it  50 

write  output  Tape  o»ioi5  >  xax(2)  .  smax 

1015  FORMAT  (  A6  ,  IHi¥  ,  F8.4  »  IH$  ) 

READ  INPUT  TAPE  0.I020*  (  AIK),  K  #  1,3  ) 

1020  FORMAT  I  3A6  ) 

WRITE  OUTPUT  TAPE  0*1025,  THMAX 
1025  FORMAT  (  20HMAXIMUM  THICKNESS  »  ,F9.6  ,  IH$  ) 

READ  INPUT  TAPE  0,1030,  I  B(K)  ,  K  A!  I  ,  5  ) 

1030  FORMAT  (  5A6  ) 

rewind  I 

CALL  SETUP!  TITLE  ) 

CALL  AXPLOT  ! XO , YO,XL ,YL ,D I VX ,D I VY , XLD * YLD *  NX ,NY , I  I ,XAX , T I TLE ! 7 ) ) 
INK  »  2 

DO  70  J  »  I  ,  JJ 
12  #  2  »  NOGEOD  +  3 

READ  TAPE  I,  S,  !  YIK)  ,  K  «  I  ,  12  ) 

SI  »  S  *  XL  /  SMAX  +  XO 
YI  #  Y!I2)  *  YL  /  THMAX  +  YO 
CALL  Plot  i  si  ,  yi  ,  i  ,  ink  ) 

70  INK  #  I 

XX  »  XO  +  XL  /  4.0 
YY  #  YO  +  YL  +  .3 

CALL  LETTER  !  XX,  YY  ,  50  *  52  ,  A  ) 

YY  «  YO  +  YL  +  .1 

CALL  Letter  i  xx,  yy,  5n  *  52  ,  b  > 

YY  #  YO  +  YL  +  2.0 

CALL  letter  !  XO  ,  yy  ,  50  ,  52  ,  TITLE  ) 

CALL  FINISH  !  30  ,  TITLE  ) 

END  FILE  8 
REWIND  I 

return 

END 
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1-^UUU  CDCQCQCDCD 
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# 


ccnplot  plot  of  final  contour  r  vs  z 

subroutine  CNPLOT  {  XO  ,  YO  ,  XSCALF  ♦  YSCALF  ) 

XSCALE  IS  SCALE  OF  X  (  Z  )  AXIS  (  !.□  FULL  SCALE  ,  .5  HALF 

YSCALE  IS  SCALE  OF  Y  (  R  )  AXIS  SCALE.  2.0  DOUBLE.  ETC.  ) 

dimension  r( 1 000 ) .z ( 1 000 ) .ak( 1 000) .F( 1 000) .X I ( 1 000) .Title! 1 2 ) . 

|RO( 100)  .ZO(  100)  .ALPHaO!  I  00) .CONS( I  00 )  .W ( I  00 )  .D(  100)  .THICK!  lOD)  . 

2RT !  I  00)  .ZT  !  I  00 )  .NC (  I  00) .DTHETa!  I  000) .FLNgThI  1000) 
dimension  a!  12)  .  Y!2ID)  .  XAX!|2)  .  YAX !  1 2 ) 

COMMON  M.R.Z.AK.F.X I  .NOGEOD.RO .ZO . ALPHAO.CONS . W . D . TH I CK  .NC . RT . ZT . 
lSMAX.RMAX.ZMAX.THMAX.JJ.TlTLE.pl .DTHETA  .FLNGTH . TSUM . FLSUM , NH I GH  . 

2  NLOW  .  DISTRT  .  ADVICE  .  SHAFT  I  .  SHAFT2 
TITLFI7)  H  602631452143 
TITLEI8)  »  602346456346 

TITLE!9)  «  645153606060 

XAX!  I  )  ff  606071  536060 
YAX! I )  #  606051 536060 

CALL  Setup  !  title  ) 

XL  »  !  ZMAX  +  THMAX  +  . I  )  *  XSCALE 

YL  if  !  RMAX  +  THMAX  +  .|  )  *  YSCALE 

DIVX  #1.0/  XSCALE 
DIVY  #1.0/  YSCALE 
NX  #  DIVX  +  .99 

NY  #  DIVY  +  .99 

call  aXPLOT  !  XO.YO.XL.YL.DIVX.DIVY. I .0.  I ,0.  NX. NY. 50. XAX. YAX  ) 
INK  #  2 

DO  280  N  #  I  .  M 
XX  #  Z!N)  *  XSCALE  +  XO 

YY  #  R!N)  #  YSCALE  +  YO 

CALL  PLOT  !  XX  .  YY  .  I  .  I NK  ) 

INK  #  I 
280  CONTINUE 
REWIND  I 

IN  I  #  2  *  NOGEOO  +  4 
IN2  #  INI  +  I 
INK  #  2 

DO  300  J  #  I  ,  JJ 

READ  Tape  I  .  S  .  !  Y!K)  ,  K  #  I  .  IN2  ) 

XX  #  Yl IN2 )  *  XSCALE  +  XO 

YY  #  Y !  IN  I  )  *  YSCALE  +  YO 

CALL  PLOT  !  XX  .  YY  .  I  .  INK  ) 

INK  #  I 

300  continue 

XX  #  XO  +  XL  /  4.0 
YY  #  YO  +  YL  +  .3 

write  OUTPUT  Tape  d.iooo.  xscale  .  yscale 

1000  FORMAT  !  9HZ  SCALE  #  .  F7.4  .  I2H  R  SCALE  #  .  F7.4  .  I  H$  ) 

READ  INPUT  TAPE  0 .  I  00 1  .  !  A!I)  .  I  #  I  ,6  ) 

1001  FORMAT  !  6A6  ) 

CALL  LETTER  !  XX  .  YY  .  50.  52  .  A  ) 

YY  #  YO  +  Yl  +  2.D 

CALL  LETTER  !  XO  .  YY  .  50  .  52  .  TITLE  ) 

CALL  FINISH  !  30  .  TITLE  ) 

END  FILE  8 
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REWIND  I 

RETURN 

END 


CPERCOV  COMPUTE  PERCENT  COVERAGE  AND  HELIX  ANGLE  AT  A  STATION 
SUBROUTINE  PERCOV  (  RR  ,  N  .  1  .  PERCNT  ,  HANGL  ) 
dimension  r( I Doo ) .z ( 1000) .aK( 1000) .f( I dod) .X I ( I  ODD) .Title!  12) . 
|R0(  I  00) .Z0(  I  00 ) .ALPHaOI I  00) .CONS ( I  DO )  . W ( I  DO  )  .D (  I  00 )  .Th I C< (  I  00  )  . 
2RT( IDO)  .ZT(  100) .Nc( IOO).DTHETa( 1000) .FLNgThI 1000) 

COMMON  M.R.Z.AK.F.XI .NOgEOD.RO ,20, aLPHAO , CONS .W . D . TH I CK .Nc . RT . ZT , 
|SMAX,RMAX,ZMAX,THMAX,JJ,TITLE.PI .DTHETA.FLNGTH.TSUM.FLSUM.NHIGH. 

2  NLOW  .  DISTRT  ,  ADVNCE  .  SHAFT  I  ,  SHAFT2 
IF  (  AK(N)  )  60  .  162  .  60 

60  FOK  H  FIN)  /  ABSF  (AKIN)  ) 

W02  »  Win  /  2.0 
RHO  n  FOK  *  RR 

IF  I  RR  -  CONSII)  )  140  .  130  .  IDO 

IDO  RMC  «  SQRTF  I  I  RR-CONSII))*  (RR+CONSII)  )  ) 

HANGL  #  ATaNF  I  CONSII)  /  RMC  )  *  180. 0  /  PI 

SINA  »  CONS! I )  /  RR 

COSA  #  RMC  /  RR 

FI  »  RHO  ♦  SINA  +  W02 

FAC  »  I  RHO  +  FI  )  ♦  (  RHO  -  FI  ) 

IF  I  FAC  )  120  .  120  .  I  10 

110  Y2  #  -  COSA  *  FI  +  SiNA  *  SORTF  I  FAC  ) 

X2  #  COSA  *  Y2  /  SINA  +  W02  /  SiNA 

DPHI2  n  ATANF  I  ABSFI  Y2  )  /  I  RHO  +  X2  )  ) 

F2  #  RHO  *  SINA  -  W02 

FAC2  »  I  RHO  +  F2  )  ♦  I  RHO  -  F2  ) 

FAC2  #  MAXIF  I  FAC2  ,  0.0  ) 

Y3  »  -COSA  *  F2  +  SINA  ♦  SORTF  I  FAC2  ) 

X3  #  COSA  *  Y3  /  SINA  -  W02  /  SINA 

DPHII  n  ATANF  I  Y3  /  I  RHO  +  X3  )  ) 

PERCNT  U  FOK  ♦  I  DPHII  +  DPHI2  )  /  PI 
GO  TO  170 

120  F3  #  RHO  *  SINA  -  W02 

FAC3  n  I  RHO  +  F3  )  ♦  (  RHO  -  F3  ) 

FAC3  »  MAXIF  I  FAC3  ,  0.0  ) 

Y3  #  -  COSA  *  F3  +  SINA  ♦  SORTF  I  FAC3  ) 

X3  n  COSA  *  Y3  /  SINA  -  W02  /  SINA 

DPHII  #  ATANF  I  Y3  /  I  RHO  +  X3  )  ) 

Y2  #  -  COSA  *  F3  -  SINA  *  SORTF  I  FAC3  ) 

X2  »  COSA  *  Y2  /  SINA  -  W02  /  SINA 
RHOX2  n  RHO  +  X2 
IF  I  RHOX2  )  122  .  122  ,  126 

122  DPHI2  n  ATANF  I  ABSFI  RHOX2  /  Y2  )  ) 

DPHI2  »  DPHI2  +  PI  /  2.0 
GO  To  128 

126  DPHI2  n  ATANF  I  ABSF(Y2)  /  RHOX2  ) 

128  PERCNT  n  FOK  ♦  I  DPHII  +  DPHI2  )  /  I  2.0  ♦  PI  ) 

GO  TO.  1 70 
130  HANGL  #  90.0 

GO  TO  160 
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140  HANGL  #  D.O 

RHOMIN  #  FOK  *  CONS(I)  -  W02 
IF  (  RHO  -  RHOMIN  )  150  «  160  »  160 

150  PERCNT  #  0.0 
GO  TO  170 

I  60  RHOT  #  FOK  *  CONS(  I  ) 

F4  #  RHOT  -  W02 

FAC4  »  (  RHO  -  F4  )  *  (  RHO  +  F4  ) 

IF  (  FAC4  )  150  .  150  >  161 

161  Yl  »  SORTF  (  FAC4  ) 

DPMI  I  ff  ATaNF  (  Y  I  /  F4  ) 

PERCNT  FOK  *  DPHI I  /  PI 
GO  TO  170 

162  IF  (  RR  -  CONS! I)  )  168  .  168  .  164 

164  RMC  »  SQRTF  (  (RR+CONS(I)  )  *  (  RR-CONS(I)  )  ) 
HANGL  #  ATaNF  (  CONS(I)  /  RMC  )  *  180.0  /  PI 
PERCNT  #  W( I )  /  (  PI  *  RMC  ) 

GO  TO  170 
168  HANGL  »  0.0 
PERCNT  #  0.0 
170  CONTINUE 
RETURN 
END 

♦label 


CADJUST  ADJUST  STARTING  HELIX  ANGLE 

SUBROUTINE  ADJUST  (  I >  NA  .  NR  »  FRACT  ,  EPS  .  LL  ) 

DIMENSION  R( I  000) > Z ( | 000 ) * AK < I  000) »F( 1000) »X I ( I  000) »TiTLE(  12 ) . 

|RO(  I  00) »ZO( I  00) .ALPHaOI I  00) »CONS( I  00 ) . W ( I  00 ) . D (  I  00 ) . Th I CK (  I  00 ) » 
2RT( 100) *ZT( 100) »NC( IOa)»DTHETA( 1000) .FLNGTHI 1000) 

COMMON  M»R,Z»AK*F»XI  »NOGEOD>RO. ZO» ALPHAO.CCNS » W . D > TH I CK  ,NC » RT ,ZT . 
|SmAX,RMAX,ZmAX,THmAX,JJ,TITlE,PI .DTHETA,FLMGTH,TSUM»FlSUM,NHIGH, 

2  NLOW  »  DISTRT  ,  ADVNCE  .  SHAFT  I  ,  SHAFT2 
CONV  #  180.0  /  PI 

iter  #  0 

AZERO  #  ALPHAOI  I  ) 

RZERO  #  RO( I ) 

ZZERO  #  ZO( I ) 

C  #  CONSI I ) 

IF  (  ALPHAO(I)  -  89.0  )  30  .  30  »  20 

20  RO( I )  #  RMAX 

ALPHAO(I)  ATaNFI  CONS(I)  /  SORTf(  R0(I)*»2  -  CONS(I)»*-2)  )  ♦CONV 
DO  22  N  #  2  *  M 

IF  (  RO(  I  )  -  R(N)  )  24  »  24  .  22 

22  continue 

24  ZO(  I )  #  Z(N) 

30  FRC#  FRACT 

RV  #  FlOATFI  NA  )  /  FLOATF  (  NB)  +  ADVNCE 
AAZERO  #  ALPHAOI  I  ) 

40  CONTINUE 

DELA#  RV  -  FRC 

IF  (  ABSFI  DELA)  -  EPS  )  110  »II0  »  50 

50  DTDA  #  0.0 

CSQ  #  CONS( I )**2 


79 


COSF  (  ALPHAO(I)  /  CONV  ) 


55 

60 


RCOS  »  RO(  I)  * 

NLI  NLOW  +  I 
NH  I  NHIGH  -  I 

SQ2  #  I.D  /  SqrTf  (  R(NL|)**2 
DTDA  #  DTDA  -  F(NLOW)  *  RCOS  * 
I F  (  NH I  -  NL I  )  85  »  55  »  55 

DO  80  N  if;  NL  I  ,  NHI 
IF  (  AK(N)  )  60  .  70  .  60 
SOI  #  SQ2 

SQRTF  (  R(N+I )**2 
+  F(N)  *  RcOS  #  ( 


-  CSQ 
SQ2  / 


) 

AK(NLOW) 


-  CSQ 

-  SQ2 


-  AAZERO 
)  -  5.0  ) 


120 


SQ2  ff  1.0/ 

DTDA  #  dTda 
GO  TO  80 

70  DTda  #  DTda  +  Rcbs  *  R(N)  *  ( 

80  continue 

85  DTDA  «  DTDA  +  F(NHIGH) 

DTDA  #  2.0  *  DTDA 
IF  (  ABSF (  DTDA  )  -  ,01  )  140 

90  DALPHA  «  DELA*  360.0  /  DTDA 
ALPHAO(I)  #  ALPHAO(I)  +  DALPHA 
IF  (  ITER  -  in  )  100  .  150  .  150 

1 00  iter  »  ITeR  +  I 

CONS(I)  #  RO(I)  ♦  SINF  (  ALPHAO(I) 
CALL  DELTHA  (  I  ) 

NLOW  #  NLOW 
NHIGH  ff  NHI6H 
FRC  #  TSUM  /  360.0 
GO  TO  40 

I  10  DALPHA  ALPHAOt I ) 

IF  (  ABSF  (  DALPHA 
I  20  LL  #  0 

GO  TO  170 

130  WRITE  OUTPUT  TAPE  6  ,  1000  »  DALPHA 
1000  format  (  20HD  CHANGE  IN  ALPHA  ,  , 

lESIC  DISTORTED  INSTEAD  ) 

GO  TO  160 

140  WRITE  OUTPUT  TAPE  6  ,  lolo  ,  dTDA 
1010  FORMAT  (  22Hn  D  THETA  /  D  ALPHA  « 
lALPHA  WOULD  BE  REQUIRED  -  GEODESIC 
GO  TO  160 
150  WRITE  OUTPUT  TAPE  6.  1020 
1020  FORMAT  (  7|H0  ALPHA  DID  NOT  CONVERGE 

IDISTORTED  instead  ) 

160  LL  fH  I 

CONS( I )  #  C 
RO( I )  »  RZERO 
ZO  (  I  )  ff  ZZERO 
ALPHAOI  I  )  ff  AZERO 
CALL  DELTHA  (  I  ) 

170  CONTINUE 

return 

END 

*LABEL 


SQI  )  /  aK(N) 


Z(N+|)  -  Z(N)  )  *  (  $02  ) 

*  RCOS  *  SQ2  /  AK(NHIGH) 

140  .  90 


/  CONV  ) 


130 


I  30 


FI  0.6  ,  43H  »  TOO  GREAT  -  GEOD 


.  F9.6  , 

DISTORTED 


7IH  »  LARGE 
INSTEAD  ) 


CHANGE  IN 


IN  10  ITERATIONS  -  GEODESIC 


CNOFACT 


alters  fraction  so  no  common  FACTORS 


80 


SUBROUTINE  NOFACT  (  NUMER  »  IDENOM  ) 
I  00  JJJ  #  0 

200  JJJ  »  JJJ  +  I 

Ml  IGCD  (  NUMER, IDENOM  ) 

IF  (  Ml  -  I  )  300  ,  300  ,  ?50 

250  GO  TO  (  I  ,  2,  3,  A,  5  )  ,  JJJ 

1  IDENOM  #  IDENOM  +  I 
GO  TO  200 

2  IDENOM  »  IDENOM  -  2 
GO  TO  200 

3  IDENOM  ff  IDENOM  +  I 
NUMER  fl  NUMER  +  I 
GO  TO  200 

4  NUMER  iff  NUMER  -  2 
GO  TO  200 

5  NUMER  ff  NUMER  +  2 
IDENOM  ff  IDENOM  +  I 
GO  To  lOD 

300  CONTINUE 

return 

END 


FUNCTION  IGCDIMM.NN) 

c  program  author  m.elson,  central  data  processing, I / I /65 

Mi«MM 
N)I(NN 

IF(M-N)2,2, I 

1  I#M 
M»N 
N#I 

2  IGCD#M 
IGCDI»XM0DF(N,M) 

IF(  IGCDI  14,4,3 

3  NSfM 
M»IGCD  I 
G0T02 

4  return 

END 

»LABEL 

CSHIFT 

subroutine  shift  (  N  ,  SHIFT2  ) 

dimension  R(  I  000  )  ,Z  (  I  OCO  )  ,  AK-(  I  000  )  ,F  (  I  000  )  ,X  I  (  I  000  )  ,TlTLE(  I  2  )  , 

I RO (  100)  ,Z0(  I  00 ) , ALPHAO (  IQO)  ,CONS(  100)  ,W(  |00) ,D(  100) ,THICK(  100) , 
2RT (  I  00  )  ,ZT (  IQO)  ,NC (  100)  ,DTHET  A(IOnn),FLNGTH(IQOO) 
dimension  R|(I00  )  ,  R2(IOa  )  ,  PHIIIOO  )  ,  XC(IOD  ) 

COMMON  M,R,Z,  Al^,F,XT  ,N0GE0D  ,  RO  ,Z0  ,  ALPHAO  ,CONS  ,  W  ,  D  ,THI  CK  ,Nc  ,  RT  ,ZT, 
I SMAX,RMAX ,ZMAX ,THMAX ,JJ ,T I TLE ,Pl ,DTHETA  ,FLNGTH ,TSUM , FLSUM , NHI GH , 

2  NUOW  ,  DISTRT  ,  ADVNCE  ,  SHAFT  I  ,  SHAFT2 
COMMON  AA,BE,CC,DEL,DELRH0,NSTART 
COMMON  LLL  ,  RHOMIN  ,  FR  ,  TMIN 
COMMON  Rl  ,  R2  ,  PHI  ,  XC 
SHIFT2  ff  0.0 
I  N  -  I 

5  IF  (  AK( I )  ) 


I  0  ,  60  ,  10 


rvonnnoonn 
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in  IF  (  PHKi)  -  9n.n  )  5n  .  5o  »  an 
an  IF  (  PHKI)  -  lan.n  )  an  »  4n  »  4n 

3n  SHiFTa  «  SHiFia  +  rki)  -  pain  *  cosf  i  phki)  *  pi  /  isn.n  ) 

GO  To  60 

40  SHIFTa  »  SHIFT2  +  RKI)  +  Rad) 

GO  TO  60 

50  SHiFTa  »  SHIFT2  +  RKI)  *  (  1.0  -  COSF  (  PHKI)  *  PI  /  180.0  )  ) 
60  IF  (  I  -  N  )  70  »  80  »  80 

70  I  »  N 
GO  TO  5 

80  continue 
return 

END 

♦LABEL 


CAXPlOT  ORAW  AXES  FOR  PLOTS 

subroutine  AXPLOT  (  X0»Y0,XL,YL,DI  VX»DI  VY,XLd,-YLd»NX.NY,  I  I  ,XAX, 

I  VAX  ) 

XO  ,  YO  IS  THE  ORIGIN 
XL  ,  YL  IS  LENGTH  OF  AXES 
DIVX.DIVY  IS  DIVISIONS  PER  INCH 
XLD.  YlD  is  length  DIVISION  REPRESENTS 

NX  >  NY  IS  DIVISIONS  TO  BE  LABELED  (  I »EVERY  DIV  ,  a»EVERY  OTHER) 
II  IS  SIZE  OF  LETTERS 
XAX  ,  YAX  IS  name  OF  AXES 

DIMENSION  XAX( |2)  »  YAX( 12) 

XXL  #  XL  +  XO 
YYL  »  YO  +  YL 

CALL  plot  (  XO  ,YYL  ,1,2) 

CALL  PLOT  (  XO  ,  YO  ,  I  ,  I ) 

IX  ff  DIVX  *  XL  +  .05 

lY  DIVY  *  YL  +  .05 

YOFFI  ff  yo  +  .04 
YOFFa  YO  -  .04 
XOFFI  »  XO  -  .04 
XOFF2  #  XO  +  .04 
DO  20  I  #  I  ,  lY 
YI  ff  I 

YYI  »  YI  /  DIVY  +  YO 
CALL  PLOT  (  XOFFI  ,  YYI  ,  |  ,  2  ) 

ao  CALL  PLOT  (  XOFFa  ,  YYI  ,  I  ,  I  ) 

CALL  PLOT  (  XXL,  YO  ,  |  ,  2 ) 

CALL  PLOT  (  XO  ,  YO  ,  |  ,  | ) 

DO  10  I  #  I  ,  IX 
XI  ff  I 

XXI  #  XI  /  DIVX  +  XO 
CALL  PLOT  (  XXI  ,  YOFFI  ,1,2) 

10  CALL  Plot  i  xxi  ,  YOFFa  ,  i  ,  i  i 

I F  (  I  I  -  5  I  )  30  ,  40  ,  50 

30  SIZE  ff  .096 
GO  TO  90 
40  SIZE  ff  .193 


82 


GO  TO  90 

50  IF  (  II  -  53  )  60  »  70  »  80 

60  SIZE  .384 
GO  To  90 
70  SIZE  #  .768 
GO  TO  90 
80  SIZE  #  1.536 
90  CONTINUE 

YOFF  #  YOFF2  -  SIZE  -  .1 
DO  100  I  #  NX  *  IX  >  NX 
XI  #  I 

XXI  #  XI  /  DIVX 
XXXI  »  XI  »  XLD 

write  OUTPUT  TAPE  O.IOOO.XXXI 
1000  FORMAT  (  F5.2  .  1H$  ) 

READ  INPUT  TAPE  0,1002,  A 
1002  format  (  A6  ) 

XXI  »  XXI  -  2.5  *  SIZE  +  XO 
IDO  CALL  letter  (  XXI  ,  YOFF  ,  I  I  ,  52  ,  A  ) 
YOFF  »  YOFF  -  SIZE  -  .1 
XX  #  XO  +  XL  /  4.0 

CALL  letter  (  XX  ,  YOFF  ,  I  I  ,  52  ,  XAX  ) 
XOFF  #  XOFFI  -  .05  -  5.0  *  SIZE 
DO  110  I  #  NY  ,  lY  ,  NY 
YI  X  I 

YYI  YI  /  DIVY  +  YO  -  SIZE  /  2.0 
YYYI  »  YI  *  YLD 
WRITE  OUTPUT  TAPE  0»I000>  YYYI 
READ  INPUT  TAPE  0,  1002  ,  A 
110  CALL  letter  (  XOFF  ,  YYI  ,  1 1  ,  52  ,  A  ) 
XOFF  XOFF  -  SIZE  -  .1 
YY  »  YO  +  Yl  /  4.0 

CALL  letter  (  XOFF  ,  YY  ,  I  I  ,  53  ,  YAX  ) 

return 

END 


MAC  I  i!^MACRO/M  100  1 

$$  M  IS  THE  number  OF  POINTS  DEFINING  THE  CONTOUR  1002 

N#0  1003 

X  I  (  I  )  AfO  1005 

'010)  N#N+I  1007 

IF(ABSF(Z(N)-Z(N+I  ) )-.0000DI  )  I D M , I  0  I  I  , I  0  I  8  1028 

ION)  IF  (R(  N+l  )-R(N  )  )  I  0  I  2  ,  I  0  I  2  ,  I  0  I  3  1030 
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1012)  K(N)»-( 10**20) 

R  I  (N)#R( N+l ) 

R2(N)«R(N) 

XC(N)#XI (N)+R2(N) 

JUMPTO/ 1014 

1013)  K(IM)»I0**20 
RI (N)»R(N) 

R2(N)«R(N+| ) 

XC(N)#XI (N)-Rl (N) 

1014)  F(N)»IO**2a 
PHI ( N) #359.9 

XI (N+l )#X| (N)+R2(N)-R| (N) 

X2(N)»XI (N+( ) 

JUMPTO/ I  050 

1018)  K(N)#(R(N+| )-R(N) ) / ( Z ( N+ I ) -Z ( N ) ) 

F(N)#SQRTF( l+K(N)**2) 

X  I (N+l  ) #X I  (N)  +  (Z ( N+l  )-Z (N) )*F(N) 
X2(N)#XI (N+l ) 

PHI (N)#ABSF (K(N)/F(N))*360 
IF( ABSF ( K(N) )-. 00000  I ) 1 040 » I  040 ♦ I  020 
1020)  IF(K(N)) 1025,1025.1030 
1025)  XC(N)#XI (N+l )-R(N+l )*F(N)/K(N) 

RI (N)#XC(N)-XI (N+l ) 

R2(N)#XC(N)-X| (N) 

JUMPTO/ I  050 

1030)  XC(N)#XI (N)-R(N)*F(N)/K(N) 

RI (N)»X|  (N)-XC(N) 

R2(N)#XI (N+l )-XC(N) 

JUMPTO/ I  050 

1040)  R2(N)#6.283I8I*R(N) 

K(N)#0 
XC(N)  #  0.0 

1050)  IF(N-M+I . I )  I  0  I  0, I 060» 1060 
1060)  N#l 
1070)  N#N+| 

IF(K(N)-K(N-I  ) )  1080,108  0,1090 
1080)  IF(N-M+I . I ) 1070, 1200, 1200 
1090)  SHIFT#0 
I#N-I 

1100)  IF(K(  I ) )  I  M  0,  M  60, I  I  I  0 
1110)  IF(PHI(I )-90 )  I  I  20,  I  120, I  I  30 
1120)  SHIFT#SHIFT+R|(I)*( |-CoSF(PHI ( I ) ) ) 
JUMPTO/ I  160 

1130)  IF(PHI(I)- I  80)11 40, 1150, 1150 
1140)  SHIFT#SHIFT+R I ( I )-R2( I )*COSF (PHI ( I ) ) 
JUMPTO/ I  I  60 

1150)  SHIFT#SHIFT+R| ( I )+R2( I ) 

1160)  IF( I-N  +  . I )  I  I  70  ,  I  I  80, I  I  80 
1170)  I#N 

JUMPTO/ I  100 
1180)  I  m-  I 
1190)  I#I+I 

X  I ( I )#X I  ( I )+SHIFT 
X2 ( I )»X2 ( I )+SHIFT 
XC( I )#XC( I )+SHIFT 


I  032 
I  034 
I  036 
I  038 
I  040 
I  042 
4  044 
I  048 
I  050 
I  052 
I  054 
I  056 
I  058 
I  060 
I  009 
I  01  I 
1013 
10  14 
1015 
1017 
1019 
I  02  I 
I  023 
1025 
110  1 
I  I  03 
I  I  05 
I  107 
I  I  09 
I  I  I  I 
1112 

1113 

1114 
1116 
1118 
I  120 
I  122 
1201 
1203 
1207 
1209 
12  11 
12  13 
1215 
1217 
1219 
122  I 
1223 
1225 
I  30  I 
1303 
I  305 
1307 
I  309 
1311 


IF( I-M+l . I  ) I  190. I  070. I  070  1313 

1200)  L l#LINE/0.0. I  0.0  1401 

N#0  1403 

1210)  N#N+I  1405 

IF(K(N) )  1220. 1240. 1230  1407 

1220)  PC(N)#P0INT/XC(N) .0  1409 

L2(N)i¥LINE/PC(N)  .ATaNGL.  (  I  80-PHI  (N)  )  14  11 

C  I  (N)»C I RCLE /CENTER. PC(N) .RADIUS.RI  (N)  1413 

C2(N)#C  I  RCLE /CENTER. PC(N) . R ADI  US , R2 ( N )  1415 

JUMPTO/1250  1417 

1230)  PC(N)#POlNT/XC(N) .0  1419 

L2(N)i!(LINE/PC(N)  .  ATaNGL.PHI  (  N)  142  1 

C  I  (N)ii!(CIRCLE/CENTER.PC  (  N  )  ,RADIUS,R  I  (  N)  142  3 

C2(N);!<CIRCLE/CENTER.PC(N)  .  R  ADI  US ,  R2  (  N )  1425 

JUMPTO/1250  1501 

1240)  L2 (N)#LINE/PARLEL .L I .YLARGE.R2 (N)  1503 

L3(N)#LINe/(POiNT/X| (N) ,0) .PeRPTO.LI  1505 

L4(N)«LINE/(POINT/X2(N).0).PeRPTO.L|  1507 

1250)  IF(N-M+I . I ) 1210. 1260. 1260  1509 

1260)  TERMAC  1511 


MAC2  #MACR0/M  2001 

$$  M  IS  THE  number  OF  POINTS  DEFINING  THE  CONTOUR  2002 

O)!(POINT/0.0  2005 

STRTWPOINT/0.  I  0  2007 

TLON  2011 

FROM/STRT  2013 

GOTO/0  2015 

DRAFT/ON  '  2009 

N#0  2017 

2010)  N#N+I  2019 

IFIPHI (N)-I80)  20  I  2 .20  I  2 .20  I  4  2025 

2012)  II»I  2027 

JUMPTO/2016  2029 

2014)  II»2  2031 

2016)  IF(K(N) )2020. 2025. 2030  202  1 

2020)  GOBACK/C2(N)  . ON . I  I . I NTOF . L2 ( N )  2023 

2022)  GORgT/L2(N)  .ON.cMN)  2101 

TLON.GORgT/C I(N).T0.II.INT0F.L|  2103 

JUMPTO/2040  2105 

2025)  DNTCUT  A2I30 

GODLTA/. I  .0.0  B2I30 

INDIRV/-1  ,0.0  C2I30 

GO/ON, L3(N)  D2I30 

CUT  E2130 

GORGT/L3(N) ,ON,L2(N)  F2I30 

2027)  GORGT/L2(N)  ,ON,L4(N)  2132 

TL0N,G0RGT/L4(N) ,TO,Ll  2134 

JUMPTO/2040  2136 

2030)  GOBACK/CIIN)  .ON , I  I , INTOF , L2 ( N )  2  107 

2035)  GORgT/L2(N)  .0N,C2(N)  2109 

TLON,GORGT/C2(N) ,TO,II , INTOF.L I  2111 

2040)  IF(N-M+I . I )2050,2070.2070  2113 

2050)  IF(ABSf(X| (N+t  )-X2(N) )-.00D00l )20I0. 2010. 2055  2115 
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2055)  GOTo/(POINT/X| (N+l ) .0) 

N»N+  I 

IF(PHI(N)-|8Q)  2Q57»2057»2059 
2057)  II#| 

JUMPTO/2060 

2059)  II#2 

2060)  IFtK(N) )2062.206A»2D66 

2062)  GOLFT/c2(N)  » ON » I  I > I NTOF ♦ L2 ( N ) 
JUMPTO/2022 

2064)  G0LFT/L3(N)  »0N,L2(N) 
JUMPTO/2027 

2066)  GOLfT/cI(N)  . ON , 1 1  , 1 NTOf » L2 ( N ) 
JUMPTo/2035 
2070)  GOTO/0 

DRAFT/OFF 

TERMAC 


2  117 
2  119 
2  I  60 
2  162 
2  164 
2  I  66 
2140 
2  I  42 
2  144 
2  I  46 
2  148 
2  150 
2  I  52 
2  125 
220  I 
2203 


Mac3 

$$ 

$$ 

$$ 

$$ 

$$ 

$$ 

3010) 

3012) 
30  15) 
30(8) 
3020) 


3030  ) 


3040) 

3050) 

3060) 

3070) 

3080) 

3090) 


3100) 
3  110) 


#MACRO/RO,aZERO,PRIMe,M,EPS 
RO  IS  THE  RADIUS  OF  THE  STARTING  STATION 
AZERO  IS  THE  HELIX  ANGLE  AT  THE  STARTING  STATION 
PRIME  IS  THE  DESIRED  NUMBER  OF  CIRCUITS  PER  PATTERN  (  A  PRIME  NO, 
M  IS  THE  NUMBER  OF  POINTS  DEFINING  THE  CONTOUR 
EPS  IS  THE  MAXIMUM  ALLOWABLE  DIFFERENCE  BETWEEN  REVOLUTIONS 
PER  CIRCUIT  obtained  AND  REVOLUTIONS  PER  CIRCUIT  DESIRED 
NSO 
N#N+I 

IF  (  R(N)-RO  )  30  I 2»30 I  8 »30  I  5 

IF  (  N-M  )  3010*3355,3355 

NA?N-  I 
J#N 

PASS# I 

DEG#l80/3. 1415927 

ALPHA#AZERO 

SINA#SINF(ALPHA) 

COSA#COSF( ALPHA) 

CONS#RO*SlNA 

IF (CONS-R ( I ) ) 3280*3040,3040 
I  #J 
I  «I-  I 

IF(R< I )-CONS) 3060*3060*3050 
IF ( CONS-R (M))3280*3 07 0*3070 
L#J 
L#L+  I 

I  F (R ( L )-cONS ) 3090 *3090, 3080 

L#L-I 

N#I 

ASEC2'#ATANF(SQRTF(  (R(N+I  )/C0NS)**2-I  )  ) 

DBETA(N)#ASEC2 

DTHETA(N)#F(N)*DRETA(N)/K(N) 

FLNGTH(N)#R2(N)*SINF(DBETA(N) ) 

N»N+  I 


IF(K(N))3I  20*31  10*3120 

DTHET A ( N ) #CONS* ( Z ( N+ I  )-Z ( N ) ) *DEG/ ( R ( N ) *SQRTF ( R ( N ) »#2-CONS*#2 )  ) 
FLNGTH ( N ) #SQRTr ( ( Z ( N+ I ) -Z ( N ) ) **2+ ( R ( N ) *DTHETA ( N ) /DEG ) »*2 ) 


300  I 
A3000 
B3000 

)  C3D00 
D3000 
E3D00 
F3D00 

3002 

3003 

3004 

3005 

3006 

3007 
30  I  I 
3012 
3010 
30  I  3 
3015 

301  7 

30  19 

302  I 
3023 
3025 

31  03 
3  I  05 
3  107 
3  109 
3  111 
3113 
3  115 
3  117 
3  119 
3  120 
3121 
3123 
320  I 
3202 
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JUMPT0/3iaD  3203 

3120)  IF(N-L)3l3n*3l4n,3l4r|  3205 

3130)  ASECI#ASEC2  3207 

ASEC2A!ATANE(  SQRTF  (  (  R  (  N+ I  ) /CONS  )  I  )  )  32  09 

DBETA(N)#ABSF(ASEC2-ASECI  )  3211 

DTHETA(N)»F(N)*DBETA(N)/ABSF(K(N)  )  32  13 

FLNGTH  (  N  )  SSQRTF  (  R  I  (  N  )*»2+R2  t  N  )»»2-2-»R  I  (  N  ) ♦R2  (  N  )  *COSF  (  DBETA  (  N  )  )  )  32  14 

JUMPTO/3ir)n  3215 

3140)  DBETA(N)#ASEC2  3217 

DTHETA(N)#F(N)*(-ASeC2)/K(N)  3219 

FLNGTH(N)i(i‘R2(N)*SINF(DBETA(N))  3220 

TSUM»0  3221 

N))!I-I  3223 

3150)  N#N+I  3225 

TSUM»TSUM+DTHETA(N)  3301 

IF(N-L)3l5n. 3160. 3160  3303 

3160)  TSUMi!(2*TSUM  3305 

RVN)i‘TSUM/360  3307 

N#0  3309 

3170)  N»N+I  3311 

IF(RVN-N)3I80. 3 190.3170  3313 

3180)  INTGER»N-I  3315 

FRACT#RVN-INTGER  3317 

JUMPTO/3200  3319 

3190)  INTGER#N  3321 

FRACTWO  3323 

3200)  N»l  3325 

PARTNtfl/PRIME  3327 

3210)  N#N+I  3401 

PARTNI^PARTN  3329 

PARTNjS^N/PRIME  3331 

IF(FRACT-PARTN  ) 3230.3260.3220  3403 

3220)  IF (N-PRIME+I  ) 32  I  0 .3250. 3250  3405 

3230)  IF (ABSF( FRACT-PARTN  ) -ABSF ( FRACT-PARTN I  )) 3250 . 3250 . 3240  3407 

3240)  N#N-I  3409 

PARTNWPARTNI  3410 

3250)  DELSPARTN  -FRACT  3411 

IF( ABSF (DEL)-EP 5)3270. 3270. 3290  3413 

3260)  DELSO  3415 

3270)  PRlNT/3. Pass. del. N.PARTN  .FRAcT. INTGER.S  3417 

RVN, ALPHA, TSUM,  CONS.  L,I  3418 

N#I-|  3430 

3272)  N)SfN+l  3432 

PRINT/3.DBETA(N)  .DTHETa(N)  .FLNGTH(N)  3434 

IF(N-L)3272. 3380. 3380  3436 

3280)  PRINT/0  3423 

Titles  minimum  radius  is  less  than  r(i)  or  r(M)  3425 

JUMPTO/338n  3501 

3290)  IF(PASS-ID)3295. 3295. 3360  3503 

3295)  CSQr¥CONS»*2  3505 

RC0S^;R0*C0SA  3507 

N#l  3508 

SQ2#I/SQRTF(R{ I+l )**2-CSQ)  3509 

SUMSO  3511 

SUM)5!SUM-F(  I  )*RC0S*Sq2/K(  I  )  35  13 
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3300)  N^S'N+I 

IF(N-L+I )33I0.33I0>3340 
3310)  IF(K(N)  )3320«3330t332Q 
3320)  SOI#SQ2 

SQ2#I/SQRTF( R(N+| )**2-CSQ) 
SUM»SUM+f(N)»RcOS*(-SQ2+SQ| )/K(N) 
JUMPTO/33nn 

333  0)  SUMi!;SUM+RcOS»R(N)*(Z(N+l  )-Z(N)  )*(SQ2*»3) 
JUMPTO/330Q 

3340)  SUM*'SUM+F(  L)^tRCOS*SQ2/K(L) 

SUM»2»SUM 

IF( ABSF ( SUM) -.0  0000  I) 3370*3370.3350 
3350)  DALPHA#DFL/SUM  *360 
ALPHA#ALPHA+DALPHA 
PASS)!(PASS+  I 
JUMPTO/3030 
3355)  PRINT/0 

Titles  ro  as  given  is  greater  than  r  max  of  part 
JUMPTO/3380 
3360)  print/0 

Titles  alpha  did  not  converge  in  ten  passes 
JUMPTO/3380 
3370)  PRINT/0 

titles  change  in  alpha  will  not  change  theta 
3380)  TERMAC 


35  I  5 

35  I  7 
3519 
352  I 
3523 
3525 
360  I 
3603 
3605 
3607 
3609 

36  I  I 
36  I  3 
36  I  5 
36  I  7 
36  I  9 
3630 
3634 
3636 
362  I 
3623 
3625 
370  1 
3703 
3705 


MAC4  #MACR0/T2ER0,J, NUMBER 

$S  TZERO  IS  THE  STARTING  VALUE  OF  THETA 

J  IS  THE  starting  SfCTiON  FOR  THE  PLOT 
$$  NUMBER  IS  THE  NUMBER  OF  CIRCUITS  TO  BE  DRAWN 

N»J 

THETA«TZER0 
ZZ#Z( J) 

RADIAN#3.  141  5927/  18a 

ZFLAG^O 

TLON 

4010)  IF(K(N) )402a*433a*405D 
4020)  KK#-I 

IF(ABSf(ZZ-Z(N) )-.00Da0l )4030*403D,4040 
4030)  RFLAG#2 
RR#R2 I N) 

JUMPTO/408n 
4040)  RFLAG#I 
RRi((R  I  (  N) 

JUMPTO/4n80 
4050)  K)<«l 

IF(ABSF(ZZ-Z(N+I  )  )-.aaOOOI  )  4070*407  0*40(0 
4060)  RFLAG^I 
RR#R I (N) 

JUMPT0/408n 
4070)  RFLAG»2 
RR*^R2(N) 

4080)  BETA»KK*K(N)*THETA/F(N) 

X0*'KK*RR*C0SF(BETA) 


4001 
A4002 
B4002 
C4002 

4003 

4008 

4002 

4004 

4005 

4006 

4007 

4009 
40  I  I 
40  I  3 
40  I  5 

40  I  7 
4019 
4021 
4023 
4025 
4101 

41  03 
4105 
41  07 
41  09 
4111 
4113 
4115 
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YO#RRJ^SINF(BETA) 

4117 

IF(KK)4090»4D90.AI  10 

4119 

4090) 

IF(N-t.+  .  1  )4I20  . 4100. 4100 

412  1 

4100) 

RE)i;RR 

4123 

BETaZ#BETA+2*DBETA(N) 

4  125 

JUMPTO/4160 

420  1 

4110) 

IF(N-I-. 1 )4I00. 41 00.4120 

4203 

4120) 

IF( RFLAG- 1.5)4130.4130.4140 

4205 

4130) 

RE#R2(N) 

4207 

JUMPTO/41 50 

4209 

4140) 

RE«R 1 (N) 

42  1  1 

4150) 

BETAZ»BETA+DBFTA(N) 

42  13 

4160  ) 

XD#KKlfRE*COSF(RETAZ) 

42  1  5 

YDiS'RE^SINFCBETAZ) 

42  17 

IF(BETAZ-PHI (N)) 4220. 4170. 4170 

42  19 

4170) 

AI)!(SINF(PHI  (N)  ) 

422  1 

Bli¥-KK*COSF(PHI  (N)  ) 

4223 

IF (ABSF(XO-XD) -.000001 )4I80. 4180. 4190 

4225 

4180) 

A2#l 

430  1 

B2#0 

4303 

D2#XO 

4305 

JUMPTO/4200 

4307 

4190) 

SLPE  #(YD-Y0)/(XD-X0) 

4309 

A2#-SLPE 

431  1 

B2»l 

43  1  3 

D2#Y0-SLPE  *X0 

43  1  5 

4200) 

DEN0M#AI»B2-A2*BI 

4317 

IF ( ABSF( DENOM)-. 00000 1 ) 4420 . 4420 » 42  1  0 

43  19 

42  10) 

XI»-B l*D2/DEN0M 

432  1 

YI#AI*D2/DEN0M 

4323 

XOREF#XO+XC(N) 

440  1 

XIREF#XI+XC(N) 

4403 

G0T0/(P0INT/X0REF.Y0) 

4405 

DRAFT/ON 

4407 

goto/ipoint/xiref.yi ) 

4409 

DRAFT/OFF 

44  1  1 

XO#KK^tSQRTF(XI**2+Yl»>2  ) 

4413 

Y0#0 

44  1  4 

BETaZ^KBETAZ-PHI  (N) 

44  1  5 

JUMPTC/4 160 

4417 

4220) 

X0REF#X0+XC(N) 

44  1  9 

XDREF#XD+XC(N) 

442  1 

GOTO/ ( POINT/XOREF  .YO) 

4423 

DRAFT/ON 

4425 

GOTO/(POiNT/XDREF.YD) 

450  1 

DRAFT/OFF 

4502 

THETAi«!KK*F(N)*BETAZ/K(N) 

4503 

IF(KK)4230. 4230.4270 

4505 

4230) 

IF(N-L+. 1 )4240. 4260. 4260 

4507 

4240  ) 

IF ( RFLAG- 1 .5)4260.4260,4250 

4509 

4250) 

ZZ#Z(N+I) 

45  1  1 

N«N+I 

45  1  3 

JUMPTO/43  1  0 

45  1  5 

4260) 

ZZ#Z(N) 

4517 

N«N-  1 

45  1  9 

JUMPTO/43 1 0 

452  1 

427D  ) 

IF ( N- I 1 ) 4300 *43 00,4280 

4523 

4280) 

IF(RFLAG-I .5)4300*4300*4290 

4525 

4290) 

ZZ)S(Z  {  N ) 

460  1 

Nfi^N-  1 

4603 

JUMPTO/43 I  0 

46  0  5 

4300) 

ZZ#Z ( N+ 1 ) 

4607 

N#N+  1 

4609 

4310) 

IF( ABSF (ZZ-Z ( J )  )-. 00000 1 ) 4320  *  432 0  *  40 1 0 

46  1  1 

4320  ) 

ZFLA6#ZFLAG+ 1 

46  1  3 

IF (ZFLAG-2*NUMBER+. 1)401 0*4010*4430 

46  15 

4330  ) 

IF(ABSF(ZZ-Z(N)  )-.DQ000l  ) 4340 *4340  *  435  0 

46  17 

4340) 

XO#X 1 ( N) 

46  19 

XD#X2 ( N) 

462  1 

ZZ#Z ( N+ 1 ) 

4623 

NNi«;N+l 

4625 

JUMPTO/436n 

470  1 

4350) 

XO*(X2(N) 

4703 

XD#X 1 (N) 

4705 

ZZ*!Z(N) 

4707 

NN»N-I 

4709 

4360) 

Y0#THETA»RA0IAN*R(N) 

47  1  1 

DY#DTHETA(N)^fR(N)*RADIAN 

47  1  3 

yd;!<yo+dy 

47  1  5 

DENOM#XD-XO 

47  1  7 

SLPE  »DY/DENOM 

47  19 

4370) 

IF(YD-R2(N) ) 4390 *4380*4380 

472  1 

4380) 

YI»R2<N) 

4723 

XI»(SlPE  *X0+YI-Y0)/SLPE 

4725 

GOTO/(POiNT/XO*YO) 

480  1 

DRAFT/ON 

4803 

goto/ ( POinT/Xi  *Yi  ) 

4805 

DRAFT/OFF 

4807 

XOffXl 

4809 

48  1  1 

YD#YD-R2 (N) 

48  1  3 

JUMPT0/437a 

48  1  5 

4390) 

GOTO/(POlNT/XO,YO) 

48  19 

DRAFT/ON 

482  1 

GOTO/(POINT/XD*YD) 

4823 

DRAFT/OFF 

4825 

THETA#YD/(R(N)*RADIAN) 

490  1 

IF(ABSF(ZZ-Z(J)  )-.OODOOl  ) 4400  *  4400  *  44 1 0 

4903 

4400) 

ZFLAG#ZFLAG+ 1 

4905 

IF (ZF LAG-2 *NUMBER+. 1)4410*44 10*4430 

4907 

4410) 

N»NN 

4909 

JUMPTO/40 1 0 

49  1  1 

4420) 

PRINT/0 

49  1  3 

Titles 

LINE  connecting  POINTS  IS  PARALLEL  TO  L2(N) 

49  1  5 

PRINT/3*N,AI  *BI  *A2*B2*D2*PHI  (N)  *X0*Y0*XD*YD*I4(N) 

49  1  7 

4430  ) 

TERMAC 

49  1  9 
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